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Abstract: In this paper we present the first analytic computation of the six-point two-loop 
amplitude of ABJM theory. We show that the two-loop amplitude consist of corrections 
proportional to two distinct local Yangian invariants which can be identified as the tree- 
and the one-loop amplitude respectively. The two-loop correction proportional to the tree- 
amplitude is identical to the one-loop BDS result of M = 4 SYM plus an additional remainder 
function, while the correction proportional to the one-loop amplitude is finite. Both the 
remainder and the finite correction are dual conformal invariant, which implies that the two- 
loop dual conformal anomaly equation for ABJM is again identical to that of one-loop N = 4 
super Yang- Mills, as was first observed at four-point. We discuss the theory on the Higgs 
branch, showing that its amplitudes are infrared finite, but equal, in the small mass limit, to 
those obtained in dimensional regularization. 
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1 Introduction 



Amidst the shadow of tremendous progress in N = 4 super Yang-Mills (SYM4) amplitudes, 
three-dimensional Chern-Simons matter (CSM) theory has recently enjoyed a quiet surge of 
interest. This reflects an interesting dual aspect of the latter: On the one hand it is a close 
cousin to SYM theory in four-dimensions and thus provides a fruitful arena to apply the 
methods that was developed there-in. On the other, while scattering amplitudes of SYM4 
theory, both perturbative and non-perturbative, are closely related to string theory scattering 
amplitudes, such relations for CSM theory are either obscure or in some cases simply absent 
as the proper correspondence is with M-theory instead. The latter is intriguing in that it 
implies that certain novel properties that is shared between the scattering amplitudes of both 
Yang-Mills and CSM may in fact have a deeper purely field theoretical origin. 

A prominent example is the M = 6 theory constructed by Aharony, Bergman, Jafferis 
and Maldacena (ABJM) [1]. Being dual to type IIA string theory in AdS± x CP 3 background, 
it is very similar to SYM4 in terms of providing an exact AdS/CFT pair. This similarity in- 
spired the discovery of many common features between the two theories such as the presence 
of a hidden Yangian symmetry [2] (or equivalently dual superconformal symmetry [3-5]) of 
the tree- and planar loop-amplitudes [6], as well as the realization that the leading singu- 
larities of both theories are encoded by the residues of a contour integral over Grassmaniann 
manifolds [9, 10]. 

In many aspects, ABJM amplitudes are simpler than its four-dimensional relative. This 
simplicity is already reflected in the fact that only even legged amplitudes are non-trivial [11]. 
Furthermore, all one-loop amplitudes consist solely of rational functions [12-15] (multiplied 
by 7r, in a natural normalization) while the two-loop amplitudes are of transcendentality-two- 
functions [16, 17]. This should be compared to transcendent ality -two- and four-functions 
for one and two- loop amplitudes respectively in SYM4. As all one- loop amplitudes can be 
conveniently expressed in terms of a basis of massive triangle integrals, whose coefficients can 
be directly computed via recursion relations [18], the one-loop amplitude for ABJM theory 
with arbitrary multiplicity is effectively "solved". 

On the other hand some properties of CSM theory, while shared with YMs theory, demand 
an alternative explanation other than the stringy origin currently available for the latter. 
Consider the color-kinematic duality [19], which leads to non-trivial amplitude relations for 
YMs and relates the amplitudes of the gauge theory to that of the corresponding gravity 
theory to all order in perturbation theory [20, 21]. For CSM, it was shown that similar 
duality, although based on three-algebra [22] , is also present for the M = 8 [23] and M < 8 [35] 
theory. While the relations implied by the duality in YMs can be traced back to monodromy 
relations of string amplitudes [24], such correspondence does not exist for CSM theory since 
the amplitudes are not directly related to any open string amplitudes in a flat back-ground. 

*At this stage it is unclear what role, if any, AdS/CFT plays in the existence of these symmetries, as explicit 
attempts at proving self-T-duality [7] on the string or supergravity side have encounter technical difficulties 
and have not been fully carried out [8] . 
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As the color-kinematic identity allows one to obtain the amplitudes of gravity-matter theory 
from that of CSM theory, 2 the fact that gravity amplitudes can be extracted from close string 
amplitudes, render the role of string theory even more mysterious. 

In this paper our main focus is the loop amplitudes of AB JM theory, in particular, the six- 
point one- and two-loop amplitude, in the planar ('t Hooft) limit. It was shown in ref [16, 17] 
that the two-loop four-point amplitude has the same functional dependence as that of the 
one-loop four-point SYM4 amplitude. This equivalence was latter shown to persist to all 
orders in e expansion [25]. As the four-point amplitude can be uniquely determined by the 
dual conformal anomaly equation [26, 27], this results states that the anomaly equation for 
both ABJM and SYM4, up to four-points, are identical. However, taking into account the fact 
that the theory is conformal, the simplicity of four-point kinematics and the transcendental 
requirement of the finite function, this result might be deemed accidental (although not for 
the all order correspondence). At six-point, it is nontrivial that the anomaly equations should 
match. Furthermore, the anomaly equations fixes the result only up to homogenous terms 
and six-point is the first place where non-trivial invariant remainder functions might appear. 
Thus the six-point computation is an important piece of data to clarify these issues. 

As ABJM theory consists of matter fields transforming under the bi-fundamental repre- 
sentation of the gauge group SU(N)fcxSU(N)_fc, the amplitude has a definite parity under the 
exchange of Chern-Simons level k o —k. More precisely and L-loop amplitude is weighted by 
a factor of (47r/k) L+1 , and hence (odd-) even- loop amplitudes are parity (even)odd. Assuming 
that parity is non- anomalous, in order for odd-loops to give an acceptable contribution it must 
compensate for its opposite parity. Since the exchanging of k — k can be translated into the 
exchange of the gauge group, this implies that a non- vanishing odd-loop all scalar-amplitude 
must pickup a minus when cyclicly shifted by one-site. This is indeed the case. 

We construct the six-point integrand using leading singularity methods. Since it was 
shown in ref. [6] that there is only one pair of leading singularity at six-points, it is straight- 
forward to construct the integrand by choosing an integral basis consisting of integrands with 
uniform leading singularities. At one-loop there are two types, the one-loop box and massive 
triangle integrals with loop momentum dependent and independent numerators respectively. 
There are two distinct combinations of the leading singularity pair, the difference and the sum. 
The former is simply the tree amplitude while the latter is the conjugate tree-amplitude, with 
even and odd sites now belonging to the conjugate multiplet, denoted as -4g r fhifted since the 
identification of multiplets are shifted by one site. We find that those two objects do appear 
in integrand. 

In ref [13-15], the one-loop amplitude was given solely in terms of triangle integrals, 
proportional to ^shifted- This is valid up to order 0(e) as the box- integrals integrate to zero 
at 0(1). Having a result at one-loop that is valid to all orders in e will be extremely important 
for the construction of the two-loop amplitude. 3 The integrated result is proportional to a 

2 Both pure Chern-Simons and gravity in three-dimensions are topological. 

3 This was already seen for the four-point amplitude [16] where the one-loop result vanishes up to 0(e), 
yet it has a nontrivial box integrand. This integrand later becomes the seed of the two-loop integrand. The 
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step function, which as we will see nicely captures the non-trivial topology of 3d massless 
kinematics. More precisely, massless kinematics in three-dimensions can be parameterized 
by points on S . For color ordered amplitudes, distinct kinematic configurations can be 
categorized by a "winding number" which can be unambiguously defined. The sign function 
then simply distinguishes the configurations with even or odd winding number, for a given 
kinematic channel. 

With the one-loop integrand in hand, we construct the two-loop amplitude by simply 
requiring that on the maximal cut of one of the sub-loops, one obtains the full one-loop 
integrand. This fixes the integrand up to possible double triangles, which are further fixed 
by soft-collinear constraints. We compute the integrals using both dimensional reduction 
regularization as well as mass regularization. This mass regulator can actually be given a 
physical interpretation in terms of moving to the Coulomb branch of the theory and giving 
the scalars a vev, similar to that used for SYM4 [28]. Interestingly, while the result for 
the individual integrals differ between the two schemes, they give, up to additive constant, 
identical results when combined into the final physical amplitude. 

Using five- dimensional embedding formalism, the tree amplitude is multiplied by five- 
dimensional parity even integrals, while the conjugate tree-amplitude is multiplied by parity 
odd integrals. Introducing the cross-ratios (only two of the these are algebraically indepen- 
dent) 

(l-3)(4-6) (2-4)(5-l) (3-5)(6-2) 



(l-4)(3-6)' 
the two-loop amplitude is given as: 
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(1.2) 



where BDSq is the one-loop MHV amplitude for M = 4 SYM [29, 30], with proper rescaling 
of the regulator to account for the fact this is at two-loops, and the remainder function Rq is 
given as 

1 



-2tt 



Liofl 



Ui) + - log Ui log u i+ i + (arccos ^fu~i) 2 



The Xi are little-group-odd cross-ratios defined in (7.4); we warn the reader that these vari- 
ables may require some care when analytically continuing to Minkowski kinematics. An al- 
ternative form of the amplitude with explicit dependence on conventional invariants is given 
in eq. (7.3). 

The presence of the BDS result demonstrates that infrared divergence and the dual 
conformal anomaly equation of the two-loop ABJM theory is identical are that of one-loop 
SYM4. Furthermore, similar to SYM4, using the mass regulator we show how the anomaly 

relevance of the 0(e) pieces can also be understand from unitarity cuts, where such terms might combine with 
collinear singularity of the tree amplitudes to give non-trivial two loop contribution. 
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equation can be converted into a statement of exact dual conformal symmetry in higher 
dimensions, with the mass playing the role of the extra dimension. 

This paper is organized as follows: In section (2) we lay out some basic conventions, while 
in section (3) we begin with the discussion of general one-loop dual conformal integrand and 
its integration in the embedding formalism. We then explicitly construct the one-loop six- 
point integrand and well as the integrated result. We end with a more detailed discussion 
of the properties of the one-loop amplitude in terms of the topological properties of three- 
dimensional kinematics. In section (4), we employ leading singularity methods and soft- 
collinear constraints to fix the two-loop integrand. In section (5) we briefly discuss two 
regularization schemes, dimensional reduction regularization and higgs mass regulation, with 
special emphasis on the latter. In section (6) we will use mass regularization to explicitly 
compute the integrals. In section (7) we combine the integrated expressions and give the 
complete six-point two-loop amplitude. We give a brief conclusion and discussion for future 
directions in section (8). 

2 Conventions 

Since we will be interested in planar amplitudes, it is useful to define the dual coordinates 

x i+ i-Xi=pi. (2.1) 

Special interest in the Xi coordinates resides in the fact that planar amplitudes in ABJM 
theories are invariant under the so-called dual conformal transformations, which act as con- 
formal transformations of the x%. To make the action of this symmetry simplest, and at 
the same time trivialize several operations which occur when doing loop computations, we 
will systematically use the so-called embedding formalism [31] (for more recent discussion, 
see [32]). 

The idea is to uplift three-dimensional Xi's to (projectively identified) null five- vectors 

yi :=(xi,l,xf) (2.2) 

such that inverse propagators become the (2,3)-signature inner product 

(i • j) ■= Vi ■ Vj ■= (xi - Xj) 2 . (2.3) 

The group of conformal transformations SO(2,3) of three-dimensional Minkowski spacetime 
is then realized linearly as the transformations of the y« which preserve this inner product. 

It was shown in ref. [6] that the tree-level amplitude and loop-level integrand in ABJM 
inverts homogeneously under dual conformal inversion: 

n 

I [An] = II \j4A n . (2.4) 
i=l 
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Due to the fact that at weak coupling the theory only has M = 6 supersymmetry, the on-shell 
states are organized into two different multiplets: 

$(77) = <f> 4 + rfipi + -e/xfi^V r <t> K + ^u K n I v J v K i , A, 

^(77) = ij 4 + rffa + -e/j^V^F + ^e/Ji^V^4, (2.5) 

where 77 are Grassmann variables in the fundamental of U(3)eSU(4). The kinematic infor- 
mation are encoded in terms of SL(2,R) spinors X a , with 

s lJ = -(ij) 2 , (y> := AfA^g (2.6) 

where Sij = xf i+2 when j = Note that x 2 j is positive when the corresponding momentum 
is spacelike, while (ij) 2 is negative in that case. For more detailed discussion of the on-shell 
variables (A", 77/ ) see ref. [2]. In this paper, we will use the convention where the barred 
multiplet sits on the odd sites. The four-point amplitude is given as [2]: 



Mim = ^ {P) ^ ) =^ {Ql) with 8\Qi):= £ niiijW,. (2.7) 



r3 &(nl 

'- with 5 2 (Q I ) := 

l<i<j'<4 

3 One-loop integrand and amplitude 



Dual conformal symmetry restricts the integral basis to be constructed of SO(2,3) invariant 
projective integrals. At one-loop, this restricts the n-point amplitude to be expanded on the 
basis of scalar triangles with appropriate numerators, as well as scalar box integrals with 
numerator constructed from the five-dimensional Levi-Cevita tensor: 

t tiiun- f 6 (<M>i>M) /o-n 
Ja [a ■ i)(a ■ j){a ■ k)(a ■ I) 

The integral in eq. (3.1) is analogous to the four-dimensional pentagon integral described in 
ref. [33] and it integrates to zero up to order e in dimension regularization [16]. To demon- 
strate how dimensional regularization is employed in the embedding formalism, we explicitly 
demonstrate this result in the following. 

We first note that eq. (3.1) can be rewritten using Feynman parametrization as 

Ibox{h j: k,l) 



- J dF e(ij,kj,dy) Jj^L (3.2) 

where dF := Y\t=i doti5{\ — J2i a i) an d Y := a\yi + a.%yj + ct-zVk + <^aVi- We now focus on the 
inner integral, which for the purpose of dimensional regularization, we define in D-dimensions: 

f 1 f d D+2 a 5(a 2 ) 1 

10 ~ m I : " r[3J J i(27r)^Vol(GL(l))(a.y)3(a./)^-3 ' (3 ' 3) 
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Let us illuminate this definition of j a by comparing it with (2.2). First, the GL(1) symmetry 
can be gauge-fixed by setting the next-to-last component of a to 1, at the price of a unit 
Jacobian. Then the 5(a 2 ) factor forces the last component of a to equal x 2 , thus reducing J a 
to the usual loop integration J ■ Finally, the factor of i is removed by the Wick rotation 
from Minkowski to Euclidean space. 

The key feature away from D = 3 is the factor (a • /) where yi := (Od, 0, 1) is the infinity 
point. This signals the breaking of dual conformal symmetry, and is required to maintain 
the projective nature of the integrand (the GL(1) invariance) for arbitrary D. This feature 
remains clearly visible when switching to the easily-obtained integrated expression: 



r [3 



(4tt)# {I ■Y) D --\\Y 2 f-^ 
Plugging this into eq. (3.2), we find that the box integral gives: 



(3.4) 



(ijkl) _[ dF r[4-f] e (i,j,M,r) r[ 3 -f] e{i,j,k,i,i) 

The first term vanishes due to the fact that Y is a linear combination of the four external 
coordinates, while the second term is at least 0(e) with D = 3 — 2e. 

As the one loop box integral vanishes, dual conformal symmetry implies that the ampli- 
tude, up to 0(e), can be solely expressed in terms of scalar triangles. However as discussed in 
the introduction, for the purpose of constructing the two-loop integrand it will be extremely 
useful (and actually essential) to have a one- loop integrand valid beyond 0(e). In the follow- 
ing, we will derive the full one-loop six-point integrand that includes both the scalar triangle 
and the tensor box integrals. We note that the form of the amplitude in terms of scalar 
triangles were given in [13, 14]. 

3.1 Leading singularity and the one-loop integrand 

At six-point there are three possible box integrals, the one mass box, two-mass-easy and 
two-mass-hard box integrals. 4 Using the five term Schouten identity of the five-dimensional 
Levi-Cevita tensor one finds the following linear identity for the box integrals: 

I box (l, 3, 4, 6) = I box (3, 4, 5, 6) + I box (4, 5, 6, 1) + I box (l, 3, 5, 6) + I box (l, 3, 4, 5) . (3.6) 

Thus the two-mass-easy integral can be expressed in terms of linear combinations of the two- 
mass-hard and one mass integrals. We will use the later two as the basis for box integrals. 
The relative coefficient of the box integrals can be easily fixed by requiring that the two 
particle cuts which factorize the amplitude into two five-point tree amplitudes, must vanish. 
Cutting in the x^-channel, shown in fig. (1), this requires four box integrals to come in the 
following combination: 

I box (3, 4, 5, 1) + I box (l, 2, 3, 4) - I box (4, 5, 6, 1) - I box (6, 1, 2, 4) . (3.7) 



Here we borrow the nomenclature of four-dimensional box integrals to denote the propagator structure. 
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Figure 1. The particular combination of tensor box integrals in eq. (3.7) combines to give vanishing 
two particle cut x 2 al = = 0. This cut must vanish as the amplitude factorizes a five-point tree 
amplitude, which vanishes. 



Using the Schouten identity, one can show that this combination is actually invariant under 
cyclic permutation by one site up to overall sign. This extra sign will be important as we will 
discuss shortly. 

The other allowed scalar integrals are the massive triangles. Their coefficients along with 
that of the boxes can be fixed by the two triple-cuts C\ t 2 (and their conjugate C 12 ), where 
the subscripts correspond to the the two distinct maximal cut, indicated as channel (1) (2) 
in fig. (2). Explicitly they are given by: 

r 3 

Ci= [J ^^^^(1,2,12,-11)^(3, 4, la,-h)M*>G>h>-h) 
^ i=i 

3 

C 2 = Jld^d^d^M-h^^M-hA^^M-hA^i) 
^ i=i 

We note that there is always an ambiguity in distinguishing C\ versus C* , since they arise from 
the two solutions of a quadratic equation. However, two convention-independent combinations 
always exist. One is the average of the two cuts C\ + C* and the other is the average of 
the leading singularities, LS\ + LS* = (C± — C*)/[4det(Zi, I2, h){Ci)\-, e.g., the numerators 
weighted by the Jacobian. Independence of the second combination follows from sign flip 
of the Jacobian on the two solutions, det(Zi, h, = — det(Zi, 1%, h)(Cl). The leading 

singularities have the following analytic form [6]: 



,nLi 



CM 



LSi = d 6 (P)S b (Q) 1 ^= 1 \ ' , LSI = LSx(+ -) . (3.8) 



2 C 25 C 41 C 63 



The functions cf, and a^ 1 are defined as 



± ._ (i\pi35\j) T ijj + 2, i - 2)(j - 2, j + 2) ±J -((-rfkih j)% ± kimnjl, m)yj) 

C^j . 2 , CM . g ) 

Pl35 Pl35 

where in the definition of cx^ 1 , the (un-barred)barred indices indicate (odd)even labels. One 
can conveniently fix the convention of C\ and C\ as: 

Ci :=2(12)(34)(56}LSi, Cf = -»■ -) . (3.9) 
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P 5 X 5 P 4 P 6 X 6 

Channel(l) Channel(2) 

Figure 2. The two maximal cuts at one-loop six-point. 



As one can check, the two combinations LS\ + LS\ and C\ + C\ both have the correct little 
group weights for an amplitude. 

The leading singularities of AB JM have a dual presentation as the residues of an integral 
over orthogonal Grassmanian [10]. As discussed in ref. [6] at n = 2/c-point there are (k — 
2)(k — 3)/2 number of integration variables in the orthogonal Grassmanian. This implies that 
at six-point, there are no integrals to be done and one has a unique leading singularity from 
the Grassmanian (plus its complex conjugate due to the orthogonal condition). This implies 
that the second maximal cut C2 and C\ must be related to C\ and C*. Indeed one can check 
that 

r, -i- r* r„ _i_ r* 

-au&ted, (3.10) 



Ci + C* C2+C2 



(12) (34) (56) (23) (45) (61) 
where we have further identified the combination as the tree amplitude rotated by one, 
Shifted (123456) := „4 tree (234561). Note that all objects in this equation have the same 
little group weights (odd under reversal of the even A's) so the identification makes sense. 

A remarkable feature of 6-point kinematics is that the expressions for C\ are explicit in 
terms of angle brackets, that is they contain no square roots. This reflects the fact that at 
six-points the cut solutions can be expressed explicitly in terms of angle brackets. Let us see 
this explicitly. At the same time, this will make apparent the following connection between 
the leading singularities and the BCFW form of the six-point tree amplitude, 

^•-Uh+LSt- ^^ -Llk + LSS, (3.11) 

in line with the original BCF logic [34] and as explained recently in [15]. The main point 
is that the on-shell condition if = Z| = in channel (1) of fig. (2) indicates that the loop 
momentum spinors can be parameterized as 

A/ 1 = Ai sin 9 + A2 cos 9, A; 2 = cos 9 — A2 sin 9) . (3-12) 

This is precisely the BCFW parameterization discussed in [6]. On the double-cut there are 
three poles as a function of cos 9, whose residues are respectively LS\, LS*, and — j4g ree . (The 
latter is located at cos# = and a computation of its residue is detailed in subsection (4.3), as 
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part of our determination of the two- loop integrand.) The desired relation then follows from 
the fact that the three residues must sum up to zero by Cauchy's theorem. For completeness, 
we record here the explicit solution which corresponds to C\ 



sm( 



^J(4e) 2 ~ (4s) 2 and cobO = -(c+)*. (3-13) 

From eq. (3.9), one also sees that Ci has a non-uniform weight under conformal inversion: 



I[Ci] 



nli J (4)444 



; I[C 2 ] 



UU J (4)444 



(3.14) 



We are now ready to use the maximal cut to completely fix the integrand. Two types of 
integrals contribute to the cut in channel (1) in fig. (2), the massive triangles as well as the 
"two- mass-hard" box integrals. As there are two solutions for the maximal cut in channel 
(1), giving different cut results C\ and C*, the massive triangle by itself cannot simultaneously 
reproduce both. This implies the need for the box integrals. On the cut the box integrals 
give: 



7^(3,4,5,1) 



2<12)<34)(56), 7^(3,4,5,1) 



2(12)(34)(56) , (3.15) 



where \c x indicates the maximal cut it is evaluated on. A simple way to verify these formulas, 
up to a common sign, is to compare their square with the square of e(a, 3,4, 5, l)/(a ■ 4) on 
the cut, using the identity 



s(ii, . . .,i 5 )e(ji, . . . , j 5 ) := det [(U -jj)], 



(3.16) 



which in fact defines our normalization of the Levi-Cevita tensor. The sign can be computed 
by a judicious use of eq. (A. 2). 

Since the one-mass box must combine with the two-mass-hard box in the combination 
given in eq. (3.7), this fixes the final integrand that reproduces the correct maximal cut to 
be (stripping a loop factor AnN/k): 



A 



A, 



tree 
6 



1-loop 



7 boa: (3, 4, 5, 1) + 7 6ra (l, 2, 3, 4) - 7 boa: (4, 5, 6, 1) - 7 6o:c (6, 1, 2, 4) 



7^(1,3,5) + 



^^7^(2,4,6). 



(3.17) 



Using eqs. (3.11) and (3.15), one can see that all maximal cut are correctly reproduced. 

An important feature of the integrand in eq. (3.17) is that it picks up a minus sign 
under a cyclic shift of the all scalar component amplitude by one-site. For the box integrals, 
this is a consequence of the linear combination dictated by the vanishing two-particle cut 
in eq. (3.7). For the triangles, this is a consequence of their coefficients: If one considers 
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the all scalar ((/> 4 (/>4(/> 4 (/>4(/> 4 </>4) component of the amplitude, from the explicit form of in 
eqs. (3.8)-(3.10), one sees that under a cyclic shift: 

C^W^^U^ = -C* 2 ^U*4>Vfa) ■ (3.18) 

These additional signs are important for a non-vanishing one-loop amplitude as we now 
discuss. In ABJM, the tree and even-loop six-point amplitudes are parity even under k — > —k, 
while odd-loops are parity odd. As parity is believed to be non- anomalous, this naively forbids 
non-trivial corrections from odd-loops unless these are odd under parity. Due to the change 
from k — > —k, we are really exchanging the two gauge group U(N)& xU(N)_fc, and thus 
resulting in a cyclic shift in the identification of the barred and unbarred- multiplet. Thus 
if the one- loop amplitude picks up a minus sign under the cyclic shift, this will compensate 
for the parity odd nature, and gives an acceptable one-loop correction. This aspect of the 
one-loop amplitude has been discussed previously in ref. [12-14]. 

3.2 The one-loop amplitude 

The box integrals integrate to zero, thus the one-loop amplitude, at order O(e ) is given solely 
by the massive triangles 5 : 



^l-loop = N ( TTjd+C*,) 7T(C 2 +C$) 

6 " k { 4^(1^7(5^3)7(1^5) 4^(2^7(4^6)7(6^2) 



The fact that (i ■ i+2) = —(ii+l) 2 motivates the following definition [14]: 

(a) 

sgn c (ij) := , ) J ' , = ±1. (3.19) 
Using eq. (3.10) the one loop six-point ABJM amplitude can thus be rewritten as 



Al~ loop =( j) -^^ifted (sgn c (12)sgn c (34)sgn c (56) + sgn c (23)sgn c (45)sgn c (61)) . 



(3.20) 



Thus the one-loop amplitude is proportional to the tree-amplitude shifted by one-site multi- 
plied by a step function. This result has been obtained previously in [12-14]. 

In closing, we note that at six-point there are only two distinct Yangian invariant, the 
sum and the difference of the leading singularity and it's conjugate. Interestingly, both 
combinations are local quantities, with the difference appearing as the tree-amplitude, while 
the sum appears as the one-loop amplitude. From eq. (3.8) this property is rather obscure, 
however due to the following non-trivial identity, equivalent to eq. (A. 4), 

+ k\l) 2 + (pi+pj+pk +pif{jk) 2 = (pi+pj +p k ) 2 (pj +Pk +Pif (3.21) 



5 The basic integral with massless internal lines, which follows easily from (3.4) with D = 3, is 

fa (a-iX.W-fc) = VISV^V^V 7 ^]- 
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one finds for example 

c 41 c 41 = ~~ 2 ' (3.22) 
^135 

Thus the denominators of the leading singularities are in fact local propagators. 6 However, 
only the sum of the leading singularities has the correct little group weights to appear in an 
amplitude. The difference does not, unless it is multiplied by sign functions, which explains 
why it can appear only at loop level. 

3.3 Analytic properties of the one-loop amplitude 

The one-loop result (3.20) displays some remarkable properties which are worth spending 
some time on. In particular, step functions are rarely seen in loop amplitudes, so we need to 
understand well why they are allowed to appear in three space-time dimensions. 

First, we would like to give some topological interpretation to the region where the 
amplitude is nonzero. In Minkowski space as null momenta can be parameterized as pi = 
£?j(l, sin 9i, cos the kinematic configuration of the scattering can be projected to a set of 
points on S . The first thing to notice is that the invariant (ij) flips sign whenever the points 
i,j on S 1 crosses each other, as was also noted in [12-14]. This is easy to see by writing the 
invariants in terms of coordinates on S 1 : 



/ Q — \ 

(ij) = 2sin( JL y A ) VEi + ie^Ej+ie. 



(3.23) 



Thus the function changes sign whenever point j crosses point i on S 1 . It is thus natural to 
divide the phase space into chambers depending on the ordering of the angles of the particles; 
the one-loop amplitude is locally constant in each of these chambers. 

Given that the product of sign functions changes sign whenever two angles cross, the 
angular dependence can be given a simple topological interpretation in terms of a "winding 
number" counting the number of angle crossings compared to the color ordering. This can be 
defined as follows: if multiples of 2tt are added to angles such that they are strictly increasing, 
< — 8i < 2tt, i = 1 . . . 7, then w := {67 — 61)/ (2vr). Then one can show 

sgn c (12)sgn c (34)sgn c (56) = k 

sgn c (23)sgn c (45)sgn c (61) 1 ' { ' ' 1 " ' 



The second factor (— l) k has a kinematical origin and originates from the factors \J [i ■ j) — ie 
which can be real or imaginary depending on whether the given channel is space-like or time- 
like, respectively. The number k then simply equals the number of positive-energy timelike 
two-particle channels. We see that the 1-loop amplitude is a highly intricate function of the 
kinematical configuration. 

As discussed in ref. [14], the fact that the one-loop amplitude is a step function can be 
readily understood from superconformal anomaly equations. Using free representation for 



The locality of the leading singularities at six-point has been recently understood as a special property of 
the orthogonal Grassmaniann [35]. 
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Figure 3. (a) The amplitude F$ in the presence of a small mass. It can be continued from < to 
8 > through a narrow window of size ~ S, which shrinks to zero size in the massless limit, (b) The 
advocated behavior in the massless setup, at a small but finite value of the coupling. A branch cut 
covers the whole imaginary axis but the discontinuity across it tends to zero at the origin. 



the OSp(6|4) superconformal generators, it was shown that acting on the one-loop six-point 
amplitude with the linear generators, one must obtain an anomalous term that is proportional 
to 5((ij)), i.e. it has support on regions where two external legs become collinear. As 
the generators are linear, single derivatives in the on-shell variables, this implies that the 
amplitude must be proportional to step functions, or equivalently, sign functions. 

However, we are rather disturbed by the notion of an amplitude vanishing in an open 
set but nonzero elsewhere — this would seem to clash with the amplitude being an analytic 
function of the external momenta. In the rest of this section, we will propose that the step 
functions behavior are not actually incompatible with analyticity of the amplitude, but are 
likely only an artifact of fixed-order perturbation theory. 

It is useful to first ask what would happen if we added small masses to the internal 
propagators still keeping the external lines massless. This could arise naturally by giving a 
vacuum expectation value to some of the scalars of the theory as discussed in section (5). In 
that case, the sign function singularity would split into two threshold singularities at = ±<5 
with 5 = y/4m 2 /EiE2- Schematically, 

sgn(0 2 - 9i) -+ F s (e 2 - X ) (3.25) 

where Fs(92 — Si) is an analytic function with an analytic window of width 25 around the 
origin. 7 This amplitude is plotted in the complex plane in fig. (3). We see that as long as 



7 The precise form of Fs can be worked out from the following exact expression for the internally massive 
loop integral, writing A = (a^-a&a;^ + m 2 (2x 2 j x 2 k + 2x%x 2 k + 2x 2 ik x 2 ik - x% - x\ k - x%)) 1/2 : 



1 _i_ . *A + m(4 + x 2 k + x% + 8m 2 ) 

„ [(a • i) + m 2 ] [(a • j) + m 2 } [(a ■ fc) + fi 2 } 8mA ° S -iA + m(x% + x 2 k + x 2 k + 8m 2 ) ' [ ' ' 

In the collinear regime x 2 j ~ m 2 , this exhibits on the first sheet a pair of logarithmic branch points at 
the threshold a; 2 , + 4m 2 = 0. However, on the second sheet there is also a square-root branch point at 
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m/0 there exists a small window of width 5 around the origin along which the amplitude 
can be rightfully continued. 

We also see clearly why such behavior is possible in three space-time dimensions but not 
in higher dimensions. In three dimensions the physical (real) phase space for a set of massless 
particles splits into chambers which are separated by singular, collinear configurations. To 
analytically continue from one chamber to the next one must avoid the singularity, since 
the amplitude is not required to be analytic around that point. But attempts to avoid the 
singularity by passing through the complex plane may fail: the singularity can be surrounded 
by cuts. 

At the massless point but at the nonperturbative level, we expect an analogous situation 
but with a nonperturbatively small window of width 5 ~ e - *^. Indeed, in a theory where 
soft and collinear quanta are copiously produced, as is ABJM, we find it unlikely for a sharp 
feature such as a sign function to remain unwashed. Rather, the backreaction of the radiation 
on the ongoing hard quanta should smear the small angle behaviour. In perturbation theory 
this would become visible through large logarithms N/klogl/6, which would have to be 
resumed at small angles. Indeed such logarithms will come out of our two- loop computation. 
Thus a more faithful model for the small angle behavior at small but finite coupling should 
be a function of the sort 

sgn((9 2 - Oi) -> ( ^ ze l (3.27) 

which can be happily continued from the left region to the right region. It would be very 
interesting to investigate the small-angle behavior quantitatively and confirm that the dis- 
continuity across the cut goes to zero as #2 — 6\ — >■ 0. 



4 The two-loop six-point integrand 

We shall now proceed to determine the two-loop six-point integrand from a variety of on-shell 
constraints. In ABJM theory we get a large number of constraints just from the fact that 
there are no 3- and 5-point on-shell amplitudes. This gives a large number of cuts on which 
the integrand must vanish. In addition, there are some very simple non-vanishing triple-cuts 
associated with soft gluon exchanges which can be used to fix the remaining freedom. 

Our first goal in this section is thus to determine the two-loop hexagon integrand using 
just the following constraints: 

0. The integrand is dual conformal invariant. 

1. Cuts isolating a five-point amplitude must vanish. 

2. Cuts isolating a three-point vertex must vanish. 

xfj — m — 'j° 2 — • The latter could be visible with physical Minkowski space kinematics, depending on 
whether the Xik and Xjk channels are time-like or not. 
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Figure 4. The triple cut of consecutive massless corners corresponds to soft exchange between the 
two external lines. In dual space, this correspond to the loop region x a approaching X{. 



3. Triple-cuts of consecutive massless corners, as shown in fig. (4), correspond to soft gluon 
exchange and must reduce to the 1-loop integrand. 

4. Absence of non-factorizable collinear divergences. 

As an example, we now show that by simply using steps 2 and 3, one completely fixes 
the four-point two- loop integrand to be that constructed in ref. [16]. This also illustrate 
the importance of obtaining the one-loop amplitude beyond 0{e). The one-loop four-point 
integrand is given by 

AT e 6(al234) 

^2 (a - l)(a-2)(a-3)(a-4)' 1 ' ' 

where the unpleasant-looking factor of y/2 is due to our normalization of the five-dimensional 
e-symbol as discussed around eq. (3.16). Now consider the triple cut of a double-box integral 
in fig. (5). On the cut, x a approaches X2 and in this limit one should recover eq. (4.1). With 
a little thought one sees that the following double-box numerator does the job: 

A ree e(al23*)(&341*) 

2 (a-l)(o-2)(o-3)(o-6)(6-3)(6.4)(6-l) 1 ' ' 

where e(a,i,j,k,*)e(b,l,m,n,*) := e(a,i,j,k, At )e(6, 1, m, n, M ). The detailed behavior of such 
numerator under the cut condition will be discussed in subsection (4.3). This however, is not 
complete as one sees that there is a non-trivial contribution to the cut (a • 3) = (a • b) = 
(b ■ 3) = 0. This separates out a three-point tree amplitude and hence must vanish. On this 
cut, using (3.16) and setting y a = to restrict to an easy subcase, the double box gives a 
nontrivial contribution 

4 tree (1 . o\2 

- ^ ±LJ1 . (4.3) 

2 (a-l)(6-l) 1 ' 

One can easily see that this contribution can be cancelled by a double triangle integral. Thus 

combining requirements (2) and (3) uniquely fixes the two-loop four-point integrand to be: 



j2-loop -^M 

•^4 — n 



f T e(al23*)(6341*) + (a-2)(6-4)(l-3) 2 



a, 6 



(a • l)(o • 2) (a • 3) (a • b)(b ■ 3)(6 • 4) (6 • 1) 



+ (sot) 



One can see that the above also satisfy requirement (1) and is equivalent to that of [16]. 
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Figure 5. The triple cut of the double box integral. As x a approaches x 2 on the cut condition, one 
should obtain the one-loop integrand given in eq. (4.1). 



4.1 Integrand basis 

We begin by constructing the most general algebraic basis of dual-conformal integrals at two 
loops. In three dimensions, the most general two-loop integral is a double-box 

T ijk;lmn {( w ,\i _ f i v l \ a )( v 2 \ &) 

2box [[V1 ' ^ ' jJ Jajb (a ■ i)(a ■ j)(a ■ k){a ■ b)(b ■ l)(b ■ m)(b ■ n) 

where v\ and Vi are some 5- vectors. Note that the presence of the numerator is required by 
dual conformal invariance, as the integrand must have scaling weight -3 with respect to both 
a and b. 8 Numerators v\ proportional to yi, yj or y^ are reducible, which would leave a-priori 
2 distinct numerators on each side. However, at 6 points constraint 2 above is very powerful 
as it requires the numerator to have zeros on any double cut isolating a massless external leg. 
For dual conformal invariant integrals, this restricts the numerators to be of the e-type 



I 



ijk;lmn 
2box 



e(a,i,j,k,*)e(b,l,m,n,*)] 



or 



ijk;kli 



2box 



e(a,i,j,k,b)] , 



where the second possibility is allowed only when (k-l) and (j -l) are both nonvanishing. Note 
that this latter parity-odd double-box integral has excessive weight on (i, k,l). At six-point, 
this can be naturally absorbed by the extra weights of C\ + C\ shown in eq. (3.14). 

Six-point double-box integrals with a three-legged massive corner, and some with two 
legged massive corners, will have two-particle cuts that factories into a product of five-point 
amplitudes as shown in fig. (6). Since five-point amplitudes vanish to all order in e, the 
contributions of these double-box integrals must cancel out on such cuts or they are not 
allowed in the integral basis. It is straight forward to see that the contributions are distinct 
and cannot cancel. Thus by imposing conditions 1 and 2 on one-loop subdiagrams the allowed 
parity even double box integrals are restricted to 



8 The absence of pentagon-boxes or more complicated topologies can be easily proved as follows. A pentagon 
would need a numerator quadratic in a. Let's assume the five external propagators involving a are oi, . . . 04 
and at- Then we can expand the numerator in terms of products (a • vi)(a • 1)2) where the (a ■ Vi) are chosen 
lie in the following basis 

(a- 1), (a -2), (a -3), {a -4) and e(al234). (4.4) 
All numerators in this basis trivially cancel some propagator, except for (e(a!234)) 2 , which would appear to 



be irreducible. However, this can be reduced using the Gram identity (3.16) together with a 



0. 
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Figure 6. Possible box integrals that have non-trivial two-particle cut which would correspond to 
factorization channel that factorize the amplitude into a product of 5-pt amplitudes. The contributions 
to the cut from each box integral are distinct, leading to the conclusion that they will not appear. 



I 2mh (i) 






e(a, i,i + 1, i + 2, *)e(6, i + 2, i — 2, i, *) 


J a,b 


(a 


i)(o • * + l)(a • i + 2) (a • 6) (6 • i + 2)(6 • i - 2) (6 • ») 


jcrab^ 


= / 




e(a, i, i + 1, i + 2, *)e(6, i — 2, i — 1, i, *) 


J a,b 


(a 


*)(o • t + l)(o • * + 2)(o • 6)(6 • % - 2)(6 • t - 1)(6 • i) 


critter 






e(a, i, i + 1, i + 2, *)e(6, i + 3, « + 4, i + 5, *) 


= / 

J a,b 


(a 


i)(o • * + l)(a • % + 2) (a • 6) (6 • i + 3)(6 • i + 4) (6 • % - 




= / 




e(a,i,i + 1, i + 2, 6) 


J a,b 


(a 


*)(o • * + l)(a • * + 2) (a • b)(b ■ i + 2)(6 • i - 2)(6 • i) 



(4.5) 

where the subscript even and odd denotes the two-mass hard integrals with parity-even and 
-odd numerators. 

The same conditions also leave box-triangle integrals 

Other choices for the numerator here, such as the other natural choice e(a,i,j,k,m), would 
be related by a Schouten identity plus double triangle integrals. The box-triangles again have 
excessive weight which would imply that they should come with factors of C\ + C* . We will 
see that they indeed arise in this way. 

Finally, the conditions applied so far leave only three double-triangle integrals, namely 

\2 



ri,i+2;i+2,i _ f + 2)' 

' / 0i6 (a-i)(o-i + 2)(o-6)(6.»)(&-* + 2) 
+ 2){i -i-2) 



2tri 

ji,i+2\i—2,i 

2tri ' Ujb ( a -i)(a-i + 2)(a-b)(b-i- 2){b-i) 

r i,i+2;i-3,i-l f (z • i + 2)(z — 1 - z — 3) 

2tn ' y ai6 (a-0(a-i + 2)(a-6)(6.i-3)(6-i-l)- K ' } 

We now finish to implement constraint 2, the vanishing of all three-point sub amplitudes. 
4.2 Constraints from vanishing three point sub amplitudes 

We consider the cut (a • 3) = (a • 6) = (b ■ 3) = which separates out a three-point amplitude 
and thus must vanish. Two types of double boxes contribute to such cut, I 2mh and I crab : and 
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they contribute: 



(1) 



(2) 



(3) 









a 


' b 





5 s(a,l,2,3,*)e(b,3,5,l,*) 



;i-3)(b-2)[(g-l)(3-5)-(a-5)(3-l)] 
(a- 2) (a - !)(&• !)(&• 5) 



i 3 i 







a 


' b 



4 s(a,l,2,3 *)e(b,3,4,5.*) 



1 ' N 5 N 

(b-2)(l-3)(a-4)(3-5) 
(o-l)(o-2)(6-5)(6-4) 



(4) 



4 


a 


b 













/ ' 3 1 



1 e(a,3,4,5,*)£(b,5,l > 3,*) 



(b-4)(3-5)[(a-5)(3-l)-(g-l)(3-5)] 
(o-4)(a-5)(6-5)(6.1) 



a 



5 e(a,l,2,3,b) 



0, 



where we've indicated the non- vanishing remainder on the cut. This was obtained by using 
eq. (3.16) to reduce the e-symbols to dot products and dropping terms which vanish on the 
cut. This can be simplified further when we take into account that the general solution to 
the cut is parametrized by y a ^ = + T a ^v, where v is any null five-vector such that v-3 = 0. 
Physically, on the cut the two loop momenta are collinear with each other. Then one finds 
that the following combination of double-box and double-triangle integrals vanish on the cut 
and are thus allowed 



L ,(i)+4« 

rl,3;3 
l 2tri 



rl,3;3,5 
L 2tri 



I 



1,3;5,1 
Itri 



/crab (1)+/ 13;3,5 



rcrittcr 



ri,i-\-2;i— 3 
1 2tri J 



T i,i+l,i+2;i-l,i-3 n T 2mh 

1 box:tri anQ 1 odd 



In addition the integral Ihnx'tri~ ' l+2 ' ? 2 is immediately ruled out. In the following, we will 
use constraint 3 to fix the relevant coefficient of the double box integrals. 
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4.3 Constraints from one-loop leading singularity 

The particular cut we will be interested in is the maximal cut of one of the sub loops with 
adjacent massless legs. This cut corresponds to a kinematic configuration where there is a 
soft-exchange between the two external legs, as can be deduced from eq. (3.12) using the fact 
that the (a • 2) only gives a pole at cos 6 — > 0. In terms of dual regions, this correspond to 
when the loop region y a —> as illustrated in fig. (4). 

More specifically, we compute the leading singularity (a • 1) = (a • 2) = (a • 3) = of 

g(a, 1,2,3,*) 
(a-l)(a-2)(a-3)(a-«)" 

Normally there are two solutions to such a cut constraint, but let us verify explicitly that 
here there is only one solution y a = y2 as claimed. To do so we expand a over a natural basis, 
such as a = a\y\ + yi + aj,yz + a\y\ + a e y e where y t := e(l, 2, 3, 4, *). Imposing the two cuts 
(a • 1) = (a • 3) = gives that a\ = and 03/04 = — (1 • 4)/(l • 3), and thus 03 oc a 2 due to the 
y 2 = constraint. Then (a • 2) ~ a 2 so the only solution is a e = 0. 

Let us work out the details. Normalizing the leading singularities in a convenient way 

F(a) 



(a ■ i)(a ■ j)(a ■ k) 

we have here 



a 



:=4 5((a.i)5((a-j))5({a-k))F(a), (4i 

residue 

d 5 a op 



4 / = 4 / d un} Q n\\ = N [ d ai da 3 da 4 da € 5{a 2 ) 
Ja J vol(GL(l)) J 

where N = \^2e(yi,y2,y3,y4,y e ) = V^y 2 . After taking the first two cuts and evaluating the 
Jacobian from the o-functions we get 

6({a-l))5((a-3)) _ N f da e 



(a- 2) ^(1-3)2(2.4)7 al 

where a(e) = y2+a t y e +a 2 2 ^4) [(Tiy^s - U 4 ] (this could be mapped to the BCFW parametriza- 
tion (3.12), since this solves the same cut constraints). Multiplying by e(a, 1, 2, 3, *)/(a • i) 
and taking the residue at a e = 0, we immediately get the box leading singularity 



e(a, 1,2,3,*) 



(a ■ l)(a ■ 2)(a ■ 3)(a ■ i) 



= V2~ 

residue 
1,2,3 



(2-*) 
(2-i) 



(4.10) 



Note that although the numerator vanishes on the cut solution y a = y2, reflecting the absence 
of three-point vertices in this theory, a nonvanishing residue remains due to the double pole 
in (4.9). The residue reflects the the exchange of a zero-momentum Chern-Simons field. 

This physical origin implies that these leading singularities are "universal", and must 
reduce to the lower-loop integrand with the loop variable y a omitted. Thus, with a normal- 
ization easily fixed from the 1-loop integrand, 



\l— loop 



= 4f~ 1Moop - (4.11) 

residue 
1,2,3 



- 19 - 



Similarly, but with an opposite sign due to k 



A l— loop 



residue 
2,3,4 



_4(*-l)-loop 



(4.12) 



These relations can easily be verified to hold for the one-loop integrand (3.17), where the 
right-hand side reduces to the tree amplitude. However, these relations must hold at any 
loop order. In a sense they are analogous to the so-called rung rule [36]. 

At two loops, this requires to see the one- loop integrand emerge on the cut, i.e. eq. (3.17). 
Indeed one finds that the various pieces of the one-loop integrand do appear from the double- 
box and the box-triangle integrands. More specifically, for the cut (a - 1) = (a -2) = (a -3) = 0, 
omitting the common v2 factor, we find the following contributions: 



jl,2,3;4,6 
box;tri 



/ 2mh m 

- 1 even \ ■*• ) 
/ Crab (l) 
jcritter^j^ 

[€(0,1,2,3,0] 



/ crab (3)^ W3,4,5,2) 



Iboxfi, 5, 1,2) 
I box (5,6, 1,2) 
J box (4,5,6,2), 
(2.t)/«(2,4,6), / 2 d m /(l) 



Itri(l, 3, 5) . 



(4.13) 



As one can see, all one-loop integrals appearing in eq. (3.17) are present. Thus the constraint 
of reproducing the one- loop integrand on the one- loop maximal cut, combined with previous 
results derived from constraint 2, fixes the two-loop integrand to be the following combination: 



^tree 



+ 



+ 



2 

Ci + Cj 

c 2 + c\ 
2V2 



i 2mh m + r 

even \ J 1 



crab 



/-Q jcritter^ 



l) + ^f' 1 



t1,3;5,1 
1 2tri 



+ cyclic 



4^(1) 



4,5,6;1,3 
boxitri 



6(0,4,5,6,1)] 



;i-s) 



+ cyclic x 2 



T 2mh 
L odd 



(2) + ff[<«. 1, 2, 3, 6)]/(2 • 6) + cyclic x 2 



+ aiI ' 



i,i+2;i— 3,i— 1 
2tri 



where cyclic x 2 implies cyclic by two sites and C\ t 2,C*2 are defined as before. The presence 
of the one-loop integrand on the cut (a • 1) = (a • 2) = (a • 3) = are shown in fig. (7). We see 
that at this point, the only remaining freedom is the triangle integrals i^t 2 '' -3 '* - . However, 
as we will now see these integrals are "badly" collinear divergent and so they are constrained 
by other physical considerations. 

4.4 Collinear divergences and the ABJM two-loop integrand 

As was demonstrated in [37] (in the context of planar J\f = 4), the exponentiation of di- 
vergences leads to constraints which can be formulated in a very simple way at the level of 
the integrand, e.g., before even performing any integral. We will now formulate similar con- 
straints in ABJM theory, but these will have a somewhat different flavor due to the absence 
of one- loop divergences. 
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Figure 7. The terms that contribute to the leading singularity (a • 1) = (a • 2) = (a • 3) = 0, which 
has to reduce to the one-loop integrand. The blue lines indicate the one-loop propagators that remain 
and uncancelled after the cut. The term in the bottom of each diagram is the numerator factor. One 
can see that the combination is precisely the one-loop answer. 



In ABJM theory, the twist-two anomalous dimensions which control the collinear and 
soft-collinear divergences begin at order (k/N) 2 , e.g. two-loops. Thus the divergences at 
two-loop are the leading ones and must be proportional to the tree amplitude in a specific 
way. Some qualitative constraints can be deduced in a simple way as follows: We note that 
soft divergences can be computed by replacing the external states by Wilson lines. Just 
this fact imposes two simple constraints. First, the coefficient of proportionality of the 1/e 2 
divergence must be a pure number, e.g. independent of the kinematics (ultimately, 6 times 
the so-called cusp anomalous dimension). Second, kinematic dependence of the subleading 
1/e divergence, which can arise from soft wide-angle radiation but not collinear radiation 
(and hence is controlled by the Wilson lines) can only be of the simple "dipole" invariants 
of the form [1 / e] log % . We will call divergences of these forms "factorizable" . These are 

^IR 

rather general constraints that any physically acceptable amplitude must possess and we will 
see that they impose nontrivial constraints on the integrand. 

We will consider the collinear divergence from the region collinear to momentum p%. 
To have a divergence we need both loop momenta to be collinear, thanks to the special 
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e-numerators, so we consider the limit 

Va^V3 + TaVA, Vb 



2/3 + ni/A- 



(4.14) 



A first requirement is that the integrals proportional to the parity-odd structure, e.g. 
Ci+C* , be finite. For the box-triangles, we find the following combination is free of divergences: 



(Ci+c?; 



4,5,6;1,3 
boxitri 



e(a,4,5,6,l)] 



(c 2 +c 2 *; 



rl,2,3;4,6 
box:tri 



e(a,l,2,3,6)] 



(1-5) J y ZJ \ (2-6) 

To see that this combination is indeed finite in the collinear region, note that in the limit 
eq. (4.14), factoring out the divergent factors l/(a • 3) (a • b)(b ■ 4) one has 

e(3,4,5,6,l) , _ e(4,l,2,3,6) 



(Ci+Ct) 



(C 2 + C* 2 ) 



(5 • l)(o • 1)(3 • 5)(6 • 6) v ZJ (2 • 6)(o • 1)(4 • 2)(6 • 6) ' 
where we've symmetrized in (a •<-)■ b). The above combination vanishes thanks to the following 
identity: 

d+Ct _ 6(6,l,2,3,4)(3-5)(5-l) 

C 2 +C* 2 e(3,4,5,6,l)(6-2)(2-4) ' 1 ' ' 

This identity is proven in appendix (A). Thus we conclude that the parity-odd part of the 
integrand in eq. (4.14) is already complete, provided that the box-triangle numerators are 
chosen as there. 

We now turn to the parity-even sector. We need to study the divergences of the yet- 
unconstrained integral l\u-f^ m more detail. Integrating out the remaining variables around 
the limit (4.14) one easily obtains the divergent contribution from the collinear region 

JL j 2 f dT a dn 

/o, 6 (a-3)(o.6)(6-4)(a.l)(6-6) K ° gM J Q<Ta<Tb<00 ^r a (n - r a )(a • 1)(6 • 6) 

where a and b are as in (4.14). Such a divergence violates factorizability in two ways: it de- 
pends on yi through (a-1) and on y§ through (b-6). This leads, for instance, to dependence on 
the cross-ratio u±. A quick look at (4.14) reveals that the only other integral with potentially 
similar dependence on T a & is J cntter . However the divergence cancels exactly, pre- integration, 
in the combination 

r rit t er (1) + / W,6 

Thus we finally arrive at the complete integrand for two-loops six-point amplitude in ABJM 
theory: 



/ 



*--(tt){ 

Ci + q 



At 



+ / crab (i) - / mttcr (i) + i£?> x - i^r - W u + cyclic 



r l,3;4,6 



+ 



+ 



2^/2 

C 2 + C\ 
2^2 



(1-5) 

1,2,3:4,1 



+ cyclic x 2 



g/(2) + ^" + cyclic x 2 



(4.16) 
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5 Interlude: Infrared regularization using the Higgs mechanism 

The two-loop amplitude is infrared divergent and must be regulated in some way. For in- 
spiration we can look at the four-dimensional sibling of ABJM, N = 4 SYM. In that theory 
there exists a canonical and self-contained infrared regularization, associated to giving small 
vacuum expectation values to the scalar fields of the theory [28]. The fields running in loops 
then acquire masses through the Higgs mechanism, rendering the loop integrations finite. 

Does a similar regularization exist in ABJM theory? As was shown in the original paper 
[1], this theory has a moduli space (C A /Zk) N where N characterizes the SU(N)x SU(N) gauge 
group and k is the level. It is described simply by diagonal vacuum expectation values (vevs) 
for the 4 scalar fields, {(f> A ) = diag(w^) (with corresponding vevs for the conjugate fields, 
(4>a) = diag((vf)^)). The identifications will play no role in what follows, although we 
will be able to see that our formulas are invariant under it. 

The first question to address is what is the spectrum of the theory at a given point 
on the moduli space. While we have not found the general answer to this question in the 
literature, this can easily be answered in perturbation theory in the usual way by studying the 
linearized action for fluctuations around the vacuum. (Due to the amount of supersymmetry, 
it is plausible that the resulting spectrum is valid for all values of the coupling, although 
this will not be important for us.) To be safe we have computed the linearized action for 
both scalar, fermion and gauge field fluctuations, and confirmed that the spectra are related 
by supersymmetry as required. These computations are reproduced in appendix (B). From 
the linearized action it is then possible to find the poles in the propagators and read off the 
spectrum. 

The result is very simple. We find that the diagonal fields remain massless, while the 
off-diagonal fields stretching between i and j acquire the mass squared 

mjj = (vi-Vi + vj-Vj) 2 - Avi-VjVj-Vi. (5.1) 

Note that this vanishes when v\ = Vj, as expected. 9 Furthermore, when mfj is nonzero, 
the computation in appendix demonstrates that the corresponding components of the gluon 
propagator lack a pole at zero momentum. Thus all modes acquire a mass. This is in contrast 
with the mass-deformed supersymmetric CSm amplitudes discussed in [11]. 

Following ref. [28], this can be used to regulate planar amplitudes. The idea is to split 
SU(N) as SU(M)x SU(N-M) with M <C N and turn on vevs only within the smaller SU(M), 
restricting attention to external states within that SU(M). Many variants are possible. For 
instance if the vev preserves the SU(M) symmetry all external states remain massless. Or a 
generic vev can break SU(M) down to U(1) M , rendering the external states massive according 
to eq. (5.1). 

However, as long as none of the SU(M) vevs vanishes, all outermost propagators in a 
Feynman diagram will be massive as depicted in fig. (8). This ensures the finiteness of the 

9 More generally, this vanishes whenever the (Zfc) iv -invariant combination Vi ® vt — Vj Vj vanishes. 
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Figure 8. Pattern of masses for the Higgsed theory following [28]. The loop propagators in the 
interior of the graph remain massless while those at the boundary, represented in bold, acquire a 
mass. External states can be chosen to remain massless or not, depending on whether the rrii are 
equal or not. 



corresponding integral. (At least for all integrals that have been considered in the literature 
so far.) 

For the purpose of regularization we will restrict to the simplest setup, taking all nonzero 
vevs to be aligned in the SU(4:) directions: vf = S^Vi. Then the mass formula reduces to 

mfj = {rrii — rrij) 2 (aligned vevs). (5-2) 

where m, := \vi\ 2 . We see that the ABJM masses behave exactly like the extra-dimensional 
coordinates in N = 4 SYM discussed in [28] ! 

In the remainder of this section, we discuss how to implement this regulator in a simple 
way within the embedding formalism. This will be applied to numerous examples in the next 
section. 

Following the extra-dimensional interpretation of the masses it is natural to enlarge the 
external five-vectors yi to six-vectors 

Vi = ( x h h x i + mi) (5.3) 

with a inner product defined such that (i ■ ip 6 ^ = for vectors of this form. Then one can 
verify that (i ■ j)^ = (xj — xj) 2 + (rrii — rrij) 2 ! , automatically generating the correct internal 
masses provided that the 6-dimensional product is used in propagators. Furthermore, the 
on-shell constraints are simply (i ■ i+1)^ = 0. Regarding loop integrations, we set the extra- 
dimensional component of loop variables to zero, a = (x, l,x 2 ,0), e.g. the loop variables 
remain 5-dimensional. Then all propagators come out correctly. 

Since only the five-dimensional components of vectors yi couple to the loop variables, it is 
immediate that all Feynman parametrization formulas in section (3) go through unchanged. 
One must simply continue to use the five- dimensional inner product (i • j) := (xj — Xj) 2 + 
m 2 + m 2 in them. The five-dimensional inner product of the external y's obeys the following 
identity 

(i ■ i+1) 2 - (i ■ i)(i+l ■ i+1) = (5.4) 
which can be seen to be equivalent to the on-shell relation (5.2) for the external states. 
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A simple consequence of this procedure is that as long as integrands written in terms of 
the j/' 6 ) are SO(2,3)-covariant, resulting amplitude will be invariant under the modified dual 
conformal generator 



K^A n = 0, where K 11 = ^ 



Xi TT-T7 — ^X^Xi-- IXiTXli— X; 



(5.5) 



This equation is essentially trivial by assumption, and will remain true as long as the integrals 
are indeed rendered finite by the regularization. 

An important question is whether the SO(2,3) symmetry is an actual property of the 
ABJM integrand even for finite values of the masses. We expect this to be the case, although 
we cannot prove it. The logic is that the dual conformal symmetry SO(2,3) is associated with 
integrability, which we do not expect to vanish into thin air just because one moves away from 
the origin of moduli space. Indeed, physically, the Higgs branch can be explored by considering 
amplitudes at the origin of moduli space but with soft scalars added, as was demonstrated 
in the context of tree amplitudes in refs. [38]. By exploring such a construction in three 
dimensions, it might even be possible to establish whether the dual conformal symmetry of 
tree amplitudes, hence presumably of loop integrands by unitarity, holds away from the origin 
of moduli space. 10 

In the present paper we work only to lowest order in the masses, e.g. we keep only 
the logarithmic dependence on them. At that level the SO(2,3) symmetry is more or less 
tautological as it is the same as the existing dual conformal symmetry. Thus to logarithmic 
accuracy in the masses the validity of (5.5) is already guaranteed by existing results. 

Let us elaborate on eq. (5.5). A consequence of it together with the on-shell condition 
(5.4) is that the dependence on the individual rrii can be determined simply from consideration 
of conformal weights. For instance, suppose an amplitude is known in the case that all internal 
masses are equal and all external masses vanish. Then, the most general "aligned" case with 
internal masses rrii (and thus generic external masses) can be obtained (to the same order in 
the small mass expansion) through the simple substitution 



Mir 



The point is that there are no ratios of the masses invariant under (5.5). Therefore, with no 
loss of generality, the Higgs regulator to logarithmic accuracy (and perhaps more generally) 
can be summarized by the simple rule 



Vi — > y% + ^irVi, e.g. (xi, l,xf) — > (xi, 1, xj + fif R ) (5.6) 



for each external region momenta, where we only need to keep track of the five-dimensional 
components of the yi, and where the massless external momenta remain undeformed. This 
recipe is the main result of this section. 



10 Dual conformal symmetry of maximal super- Yang-Mills at finite values of the masses can also be established 
by considering the symmetry as a property of the higher dimensional parent theory [39, 40]. 
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Alternatively, infrared divergences can be regulated using dimensional regularization. 
Due to the presence of Levi-Cevita tensors in Chern-Simons theory, dimensional regulariza- 
tion has always been used with a great deal of caution. However, as one can use tensor 
algebra in three-dimensions to covert the Levi-Cevita tensors into Lorentz invariant scalar 
dot products and analytically continue to D = 3 — 2e, this is more similar to dimensional 
reduction regularization, commonly applied to supersymmetric theories. This regularization 
scheme has been shown to be gauge invariant up to three-loops for Chern-Simons like theories 
in ref. [41], and has also been applied to Wilson- loop computations in ref. [42, 43] establishing 
duality with the amplitude result. We will demonstrate below that the individual dimension- 
ally regulated integrals differ from the mass regulated result functionally. However, when 
combined into the physical amplitude, the two regulated results agree. 

6 Computation of two-loop integrals 

In this technical section we describe our computation of the two-loop integrals relevant for 
the two-loop hexagon. Although we feel that some of the tricks employed here can find 
application elsewhere, the reader not interested in these details can safely skip to the next 
section. 

Certain integrals, or combinations of integrals, are absolutely convergent and can be com- 
puted directly in D = 3. Examples are I%&t, jcntter^-g _|_ J^ 4 ' 6 , or a certain combination of 
the odd box-triangles described below. It is very convenient to treat these combinations sep- 
arately since they can be evaluated without regularization. Furthermore they automatically 
give rise to functions of the cross-ratios U{. On the other hands, IR divergent integrals cannot 
be avoided. In this section we will use the Higgs regulator described in the previous section. 

We will describe the various steps in our integration method, starting from the steps 
common to all integrals. 

6.1 A Feynman parametrization trick 

There is a particular version of Feynman parameterization which is particularly effective for 
our calculations. It is inspired by a formula obtained in [44] using intuition from Mellin space 
techniques, but can be derived very simply by a judicious change of variable in standard 
Feynman parameter space as was demonstrated in ref. [45]. 

1 3 - 3 1 

We illustrate it in detail in the case of the double triangle I 2 l r \ ' , the other integrals 
are entirely similar. The first step starting from the definition (4.7) is the usual Feynman 
(Schwinger) trick 



where A = Ya=i,3 a iVi and B = £i=i,3 b iVi- 

A word about the notation. The l/vol(GL(l)) symbol means to break the projective 
invariance (oj,6i) —> a(ai,b{) by inserting any factor which integrates to 1 on GL(1) orbits. 
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Standard choices include 5(^2 a>i — 1), which give the Feynman parameter measure dF in 
section (3), or 5{a\ — 1) which give rise to Schwinger parameters. The choice of a gauge-fixing 
function will play no role in what follows, and for all practical purposes it can be gleefully 
ignored it until the final step. 

The loop integrals over a, b can all be done using only the one- loop integral (3.4). Since 
the Higgs regulator already renders the integrals finite, we set D = 3 immediately to obtain 
(to avoid cluttering the formulas in this section, we will strip a factor for each loop) 



By repeatedly using this formula and its corollary valid for b 



2 



i r°° r i 11 

r (3) / df / — — — — = --== , (6.3) 



a {a-AY{a-b) w 7 J J a (a-{A + fb)f 2 fi^A{A-b) 



we derive the following, key formula: 

1 If 1 



(a • Af{a ■ b){b ■ 2 J b . b)(b . Bf 

de 



oc 



\A- A{\{B + eA) ■ (B + eA)) 3 / 2 

dc f'°° 1 

de ,, t , , i, ' , , ; „^, -(6-4) 



torVcJo (c\A-A+\{eA + B)-{eA + B)) 

The key idea here is the introduction of the new Feynman parameter c in the last step, as 
done in [45]. Although it could be removed immediately, it will prove advantageous to leave 
it untouched until the final stage. For example, this will allow us to postpone dealing with 
square roots until the very end. 

Upon substituting (6.4) into (6.1), one notes that the variable e is charged under both 
GL(1) symmetries. Therefore, it is allowed to gauge-fix one of them by setting e = 1, which 
effectively locks the two GL(1) together. This will always be the case: the variable e is always 
removable in this way. Thus we have 

,13531 = /[^OlOsMs] (1-3) 2 

2t " J An^-cJ vol(GL(l)) ( {1 + c y. A ■ A + A • B + ±B ■ B) 2 ' 1 ' J 

As mentioned, this is similar to the formula for the double-box obtained in [44]. 

So far all we have done is rewrite the standard Feynman parameter integral in some 
specific form. As we will now see, in all cases the variables at, bi can be integrated out rather 
straightforwardly, and will generate some logarithms or dilogarithms to be integrated over c. 
The c integration in the final step then poses no particular difficulty. 
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6.2 Divergent double-triangles 

Let us see carry out the remainder of this procedure for I^trf 1 starting from (6.5). 

Notice that we haven't said anything about the regularization yet. This is because ev- 
erything is fully accounted for by the rules (5.6). According to it, we simply have to take 
(i-j) ^4 + 2/4,. 

After evaluating the dot products and doing a simple rescaling of the integration variables, 
the double-triangle is thus easily seen to depend only on the ratio e := 



ri3;3l = f°° dc I" [d 3 aia 3 bib3] 
2t " J AttVcJ vol(GL(l)) 



Att^c J vol(GL(l)) [(ai+&i)(a 3 +&3) + ca 1( i 3 + e ((ai+a 3 +6i+6 3 ) 2 + c( ai +a 3 ) 2 ) 

We are interested in the small mass limit e <C 1. The only sensitivity to e comes from the two 
regions where a*i, b\ ~ e or 03, 63 ~ e; these correspond physically to collinear configurations. 
However, everywhere else we can ignore e. Consequently, let us parametrize the variables as 

1 h h [d 3 aia 3 bib 3 ] 

ai = 1, 61 = x, a 3 = a, 6 3 = ay, — = adadxdy 

vol(GL(l)J 

such that the dangerous regions are a = and a = 00. Since the region a > 1 contributes 
the same as a < 1, by symmetry, we need only consider the former and multiply it by 2. 
Furthermore, in that region, we can neglect a in the terms proportional to e since they are 
only needed when a — )■ 0. Thus 



4rf = 2 / — - / ada 



dc f 1 , [°° dxdy 



4ttvW J (a((l + x)(l + 2/) + c) + e((l + x) 2 + c)) 2 
1np . (l+s)(l+y)+c 1 2 
' otedy- - - - = ~ lo S + ClWRj- (6-6) 



WcVo ((l + x)(l + y) + c)' ''!:; 



(This was most readily done by evaluating the integrals in the following order: y, x and c.) 

The second type of double-triangle I 2 & ' is entirely similar. The general formula (6.4) 
then gives directly, after a simple rescaling of the variables, 



r l3;35 f°° dc f [d 3 ai a 3 b 3 b 5 } 
2tri - J Any/c-J vol(GL(l)) 



X 



{(ai+b 5 )(a 3 +b 3 ) + ai& 5 + ca x a 3 + e' ((ai+a 3 +6 3 +6 5 ) 2 + c(ai+a 3 ) 2 ) ) 2 



with e' = ^^l 5 . The dangerous region is the collinear region a\ — > and 65 — )■ 0, so up to 

X 13 Z 35 

power corrections in e' we can drop at and 65 in the terms multiplying e'. These can then be 
easily integrated out, leaving: 

7 ^00 lr.tr (l+t>3)(l+63+c) 1 ,22 

r 13;35 / ^ ^g (1+b3) 2 +c " fog € _ l 
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6.3 A dual conformal integral: / critter 

As our next example, we turn to the integral / cntter (i). This integral is collinear divergent, 
but it becomes absolutely convergent after combining it with I^f' 6 as explained in section 
(4). Therefore, we will only consider the sum 

ter (1) ;= / cnt t e r(1)+/ l,3;4,6 

[d 2 a ia2 a 3 } [d 2 6 4 6 5 6 6 ] f e(a, 1, 2, 3, *)e(6, 4, 5, 6, *) + (1 • 3) (4 ■ 6) (a • 2)(6 • 5) 



± 2tri 

oo r^2„ „ „ l \j2r 

iiii-ii,! i i lit i ■' . t ;-k i i i i > i .1 ! > i i 1 - . i i i i ■ r i i i i i • ' i i ) j ■ . t i 

o vol(GL(l))vol(GL(l))y aj6 (a-A)%a-b)(b-Bf ' (6 ' 8) 

Note that we have combined the two integrals into a common Feynman parameter integral, by 
inserting the inverse propagators (a. 2) (6.5) into the numerator of the double-triangle. This 
allows us to immediately set the regulating masses to zero, since we are dealing with an 
absolutely convergent integral. 

To apply the formula (6.4), we use the familiar fact that numerators turn into derivatives 
in Feynman parameter space. Thus for instance for the first term in (6.8) 

iw 7X7Z D ^3 = e ( d ^ 1, 2, 3, *)e[d B , 4, 5, 6, *) 



(a- A) 3 (a-b)(b- B) 3 v A ' ' ' ' ' x D ' ' ' ' ' (a ■ A) 2 (a ■ b)(b ■ B) 2 ' 

Thus, after setting e = 1 in (6.4) to remove one of the GL(1) symmetries as done previously, 
we obtain 

/ c " tter (D= fA f ^^^^,1,2,3,^,4,5,6,*) 



4vr^ Jo vol(GL(l)) 

+ (1 ■ 3)(4 ■ 6)(2 • d A )(5 ■ d B )) ~ -7— — — a : ^ • (6-9) 

/ {{c + l)\A- A + A-B + \B-B) 

To proceed from here, we simply integrate over the variables aj,6j one at a time. This 
can be done in an essentially automated way using the method described in detail in a four- 
dimensional context in [45]. The idea is that at each stage the integral can be decomposed 
into a rational factor which takes the form dx / (x — Xi) n with n > 1, times logarithms or poly- 
logarithms with arguments that are rational functions of x. Such integrals can be performed, 
at the level of the symbol, in a completely automated way. After this is done, we integrate 
the symbol and obtain the c-integrand as described in [45]. 

We applied this method, doing the integrals in the order 02,65,01,66 and 03, to obtain 
the symbol of a function to be integrated over c. After a step of integration by parts in c 
to remove degree-three components, we obtained the symbol of a degree-two function, which 
could easily promoted to a function 



rcritter 



1 



dc ^ — Li 2 (l - u\(c + 1)) - Li 2 (l - u 2 ) - Li 2 (l - u 3 ) - logn 2 logu 3 



o Att^/c c + 1 

^Li 2 (l - u 2 ) - ^Li 2 (l - U3) - ^logu 2 logii3 - (arccos Vul) 2 + y. (6.10) 
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Here all non-constant terms come out of the symbol computation, while the 7r 2 /3 term is 
a beyond-the-symbol ambiguity. We have fixed it by an analytic computation at the sym- 
metrical point u\ = U2 = U3 = 1, where the integral simplifies dramatically. Assuming the 
principle of maximal transcendentality for this integral, this is the only possible ambiguity. 
As a cross-check, we have verified that this result agrees with a direct numerical evaluation of 
eq. (6.9), to 6 digit numerical accuracy at several random kinematical points with Euclidean 
kinematics, which we take to confirm our assumptions. 11 



r2mh 

L even 



6.4 Another divergent integral: l£ 

We now consider a somewhat more nontrivial divergent integral, 

f2mh m [°° dc f [cfa^ashhh] (e(d A , 1, 2, 3, *)e{d B , 3, 5, 1, *)) 

y o 47tV - c j vo1(gl(1)) (( c+1) i A . A + A . B+ i B . B f lD - j 

where the derivative operators are understood to act on the rational function underneath, to 
avoid an unnecessary lengthening of the formula. 

This integral requires regularization, and as in the rest of this section we use the Higgs 
regularization described in section (5). The procedure has the following precise meaning here. 
In both the numerator and denominator, we use the shifted five- vectors defined in eq. (5.6), 
so the formula amounts to (i ■ j) — > (xi — Xj) 2 + 2/z 2 R . 

A first observation is that in all divergent regions 65 — > 0. Thus we can drop 65 from terms 
multiplying the mass in the denominator, which allows us to integrate out 65 explicitly: 12 



/ 2mh m 



dc f [d 4 aia 2 a 3 &i63] [02(2 • 5) - 2((A + B) • 5)] (2 • 5)/(l • 3) /((A + B) • 5) ; 



AnV~c J vol(GL(l)) ((1 + c)aifl3 + ai63 + ^ + ^ + ^_ x y 



where X := (£ a + ^ b) 2 + c(£ a) 2 . 

To proceed further, we need a small bit of physical intuition about this integral. It has 
collinear divergences in the region 03, 63 — > (both loop momenta collinear to p\) and in the 
region a%,b\ — > (both loop momenta collinear to p2)- In addition, there are soft-collinear 
divergences where these two regions meet. Thus a reasonable strategy is to subtract something 
which has the same divergent behavior as /i 2 R — > but which is simpler to integrate. A good 
candidate is 

:= I?Z*(1){X — ► 4(1 + c)) (6.12) 



11 Numerics with this level of accuracy can be easily obtained starting directly from (6.9) and performing the 
d2, 65 and c integrals analytically, which are readily done using computer algebra software such as Mathematica. 
The remaining 3-fold numerical integration poses no particular problem. 

12 Strictly speaking the numerator derived from eq. (6.11) contains terms proportional to /if R . However, 
due to the special properties of the e-symbol numerators, one can see that these terms only give rise to 
power-suppressed contributions. That is, they are never accompanied by compensating 1/pni power infrared 
divergences which would render them relevant. We have verified that the same is true also for the integral 
7 crab considered below. 
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since this remains finite and has identical soft and soft-collinear regions. But thanks to the 
simplified denominator, this can be integrated more easily. Indeed after a shift a\ + b\ — > b±, 
0.3 + ^3 - > b 3 together with a simple rescaling of the variables, it can be seen to depend only 

x^ x*^ 

on a single parameter e := — 1 " 1 | 35 : 

r2mh/ n x _ [°° dc f [rf 4 aia 2 a 3 fei63] a 2 + 2bi + 2b 3 

evenU ~ J 47r v Wa 1 <6i vol(GL(l)) (a 2 + b 1 + b 3 ) 2 (b 1 b 3 +a 1 a 3 c + ea 2 (l + c)) 2 

a,3<03 

log 2 e + 0(e). (6.13) 



a-i<b 3 

7vr 2 1. 2 
lnr 

12 4 



From this point we omit further details on the computation of integrals, as they proceed using 
the same strategy as in previous examples. It remains to correct for the error introduced by 
eq. (6.12) in the hard collinear regions. At fixed a 2 , ~ 1 one can see that the region 

a 3> ^3 — * produces a logarithm whose cutoff depends on X. The error is given by the change 
in the logarithmic cutoff 



7T 2 



4vr^/c J ai<bl vol(GL(l)) 61(61 + a lC ){a 2 y + &i) 2 
Li 2 (l - y) (6.14) 



6 

so that 

= + ^O^/^is) + / coll (x 2 5 /x 2 5 ). (6.15) 

This gives the result quoted in appendix (C). 

6.5 The integral 7 crab 

The final divergent integral we have to compute is 

I crab m = f°° dc [ [ rf5a i Q 2Q3&5&6fri] (e(d A , 1, 2, 3, *)e(d B , 5, 6, 1, *)) 

U y 4vr^y vol(GL(l)) (( c + i) | A • A + A • B + \B ■ B) 2 ' K ' ' 

Its evaluation is extremely similar to that in the previous subsection. The regions which 
diverge as /i 2 R — > are the |?6-collinear and pi-collinear regions, and their intersection, the 
soft-collinear region A, B — > y%. Therefore, if we denote the // 2 R -containing terms in the 
denominator by ^\^X, we see that we can neglect a 3 and 65 in X: 

X = (ai + a 2 + b 6 + bi) 2 + c{a x + a 2 ) 2 . 

Then we proceed as in the previous example: we replace the integral by the simpler one 

/C rab/(-q . = jcrab^)^ ( Gl + b X ) 2 + Cdf) , (6.17) 

which we have been able to evaluate as (setting e\ := 4/ ^ IR ^ 35 \ 

/C rab/ (1) = _ x _ g + l (1 + logU 3 )logei - llog 2 ei + lLi 2 (l - l/u 3 ). 
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The error introduced by (6.17) is by construction localized to the hard collinear regions 
and turns out to be given by the same eq. (6.14): I crab (l) = J crab/ (l) + I coll (xj 5 / 'x| 5 ) + 

13/ J '36-' 



J (xfg/^lg). Collecting the terms gives the result recorded in appendix (C). 



6.6 Parity odd box-triangles 

As shown in section (4), parity odd box-triangles appear in the six-point amplitude only in 
absolutely-convergent combinations of the form 



(1) :=(l-4)(3-6) 



(2-4)^' 6 [e(a,l,2,3,6)] (3 • 5)igS> 3 [e(a, 4, 5, 6, 1) 



+ 



e(l,2,3,4,6) ' e(3,4,5,6,l) 

The Feynman parametrization formula (6.4) reads in this case 



. (6.18) 



odd 
box\tri 



(1) 



dc f [d 5 aia 2 a 3 &4&5&6] /e(9^, 1, 2, 3, 6) 



o 47r \/c 



vol(GL(l)) 
+ (2-S A )(3-5) 



e(l,2,3,4,6) 
e(0 fl ,4,5,6,l) 



(2-4)(5-9 B ) 



(l-4)(3-6) 



e(4,5,6,l,3) 7 (( c + 1 )i A 2 + A-B + lB 2 



.2" 



where we have combined the two integrals under a common Feynman parameter integral sign. 
This may now be evaluated directly without regularization. One can see that the integrations 
over a2, 65 do not produce any transcendental functions, which suggests to do them first. In 
the process, the e-symbols neatly cancel out: 



rodd 

box;tri 



1) 



dc f [^01036466] 



a 3 (3 • 5) 



64(2 -4) 



Airy/cJ vol(GL(l)) Vai(l-5) + a 3 (3-5) 6 4 (2 • 4) + b 6 (2 • 6) 

(l-4)(3-6) 



>2 ' 



((c + l)oia 3 (l • 3) + 0164(1 • 4) + a 3 6 6 (3 • 6) + 6 4 6 6 (4 • 6))' 
The three integrations over a«, 6j remain elementary, and we obtain a pleasingly simple result 



7-odd (i\ 
- L box;tri\ ± ) 



dc log ^ log(«i(c+l)) _ log ^ arccos (y^I) 



(6.19) 



4vrVc l-ni(c+l) 2^/tti(l -iti) 

7 The six-point two-loop amplitude amplitude of ABJM 

We now construct the final integrated result. We first consider the parity even part, i.e. terms 
in eq. (4.16) proportional to . We begin by summing all the divergent integrals, or more 
specifically, ELiC^C*) + / crab (i) + I^'^ ~ ^ir?'^' 2 ), usin g the formulas recorded 
in appendix (C). This gives 



1 



log 2 



E 

i=i 



BDS 6 - vr . 



log 2 



V+3^IR 



i+U+3 



C l,i+2 X 'i+l,i+3 



t=l 



,r 



i-l,i+2 



-Li 2 (l - u») + log log 



(7.1) 
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Pleasingly, we find that all terms of non-uniform transcendentality have canceled in the sum! 
Furthermore, the sum gives nothing but the BDS Ansatz [29, 30] in the Higgs regulator!. That 
is, up to the constant term and the substitution — > 4/ij R ! The evaluation of the parity 
even terms will be complete upon adding the dual conformal integrals — ^^ = i(^ critter (*) + 



i,i+2;i+3,i— 1\ 
2tri 



As a cross-check on our evaluation of the integrals, we have evaluated the above combina- 
tion of integrals using dimensional regularization, which has been successfully implemented in 
obtaining the two-loop four-point result [16, 17]. While we find that the results for individual 
integrals differ functionally, we find perfect agreement for the combination just considered, 
up to an expected scheme-dependent constant. This constant is given in (D.3). 

We next consider the parity odd part, i.e. terms in eq. (4.16) proportional to the sum of 
one-loop maximal cuts. The I^^ 1 integrals integrate to zero at order 0(e), thus we have: 



Ci + q/S;'' 3 [^4,5,6,l)] C2 + C */™ 6 [e(a,l,2,3,4)] 



2\/2 



(1-5) 



+ 



2^/2 



(2-4) 



+ cyclic x 2. 



Using the identity (A.l) together with the definition (6.18), this can be expressed in terms of 
the dual conformal invariant finite integral 



eh + q 
2V2 

Ct + Cf 



e(3,4,5,6,l) 



+ cyclic x 2 



2V2 



;i-4)(3-6)(l-5)(3-5) 

e(3,4,5,6,l) log ^| arccos(/^ 



(l-4)(3-6)(l-5)(3-5) 2 ^M1 



til) 



+ cyclic x 2 



(7.2) 



Due to Yangian invariance, the coefficients of the transcendental functions must be expressible 
in terms of the leading singularities LSi and LSI defined in section (3.1). This can be verified 
thanks to the remarkable identity (A. 5), together with (3.10): 



(C 1 +C 1 *)e(3,4,5,6,l) 



2^/2(1 • 4)(3 • 6)(1 • 5)(3 • 5)^(1 - Ul ) 



^shiftedSgn c (12)sgn c (45) 



((34>(46> + (35>(56» 



'«34)<46> + <35><56))< 
Thus combining everything, we find the two-loop six-point amplitude of ABJM to be 



A 



2-loop 
6 



A, 



tree 



BDS 6 + Re 



+ 



^tree 
^.shifted 



sgn c ((12))sgn c ((45)) — p = I = i = 0= log — arccos( A /ui) + cyclic x 2 | 



((34)(46} + (35}(56))' 



113 



(7.3) 



where the "remainder" function Rq is given as 

3 



R 



6 



i=l 



Li 2 (l - ui) + ^loguj \ogu i+ i + (arccos yjui) 2 
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We like to stress that the remainder function, up to an additive constant, is given entirely 
by the dual-conformal finite integral J critter (<j) for the parity even structure (and by Ihox-tri^) 
for the parity-odd structure). This is in contrast with M = 4 SYM, where the remainder 
function is mixed with BDS and spread across a number of divergent integrals. 

The part proportional to .^shifted can ^ e written in a more compact way if we assume 
certain restrictions on the kinematics. We will assume so-called Euclidean kinematics, e.g. 
all non-vanishing invariants (i ■ j) are spacelike. (This is a nonempty region even for real 
Minkowski momenta.) In that case, it is correct to naively rewrite the original expression in 
terms of angle-brackets: 



arccos( A /ttT) 1 g {^-i^i^ J (1 • 4) (3 • 6) log xi 



^/ux{l-ux) ' 2i ^ Ul (l- Ul ) 2i(12><45>«34><46> + <35><46» 

where 

= (12)(45)+i((34)(46) + (35)(56)) 
Xl " (12)(45)-i((34)(46) + (35)(56))' 1 ' 

Indeed, in that region, u\ > and the first expression is always real and positive. This is also 
the case for the second expression, as can be seen from the fact that (12) and (45) are real, 
while ((34) (46) + (35)(56)) is either real or smaller in magnitude than (12)(45). Note that, 
defining cross-ratios \2 {X'i) from cyclic shifts by minus 2 (plus 2) of this expression, it can 
be shown that X1X2X3 = 1- 

This allows us, in these kinematics, to simplify the answer to 



^2-loop _ ( fV \ j ^ 



BDS e + R 6 



6, shifted 

4i 



u 2 

log — log xi + cyclic x 2 

"3 



(7.5) 



While strictly derived from eq. (7.3) in Euclidean kinematics, we expect this expression to 
be valid in other kinematic regions for a suitable analytic continuation of the variables Xi- 
Note that as Xi ~^ X^ 1 under the Z2 little group transformation of any external leg, the little 
group weight of log xi is exactly what is needed to compensate the little group mismatch of 

^tree 

^.shifted' 

8 Conclusions 

In this paper, we construct the two-loop six-point amplitude of ABJM theory. The result can 
be separated into a two- loop correction proportional to the tree amplitude, and a correction 
proportional to the shifted tree amplitude, which are distinct Yangian invariants. The correc- 
tion proportional to the first is infrared divergent and we use mass regularization. The result 
shows that the infrared divergence is identical to that of J\f = 4 super Yang-Mills and is thus 
completely captured by the BDS result. This establishes that the dual conformal anomaly 
equation is identical between the one-loop SYM4 and two-loop ABJM, which was first ob- 
served at four-points and we conjecture will persists to all points. The correction multiplying 
the shifted tree amplitude is completely finite. 
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As a comparison, we also computed the divergent integrals using dimensional reduction. 
We find that the individual integrals give different functional answer between the two regu- 
larization schemes. However, when combined into the amplitude, they give the same result 
up to a physically expected constant. 

We find in addition to the BDS result a nonzero (dual-conformal invariant) remainder 
function. This implies that the six-point ABJM amplitude cannot be dual to a bosonic 
Wilson- loop, which only captures the BDS part [43] 13 . This does not rule out a possible 
duality with a suitable supersymmetric Wilson loop, however. The reason is that, if SYM4 is 
to be of any guidance [46], the correct Wilson loop dual for amplitudes with n > 6 particles 
should reproduce, at lowest order in the coupling, the n-point tree amplitude 14 . Since no 
candidate Wilson loop with this property, or even just the correct quantum numbers, are 
presently available in the literature, we find it hard to say anything conclusive about the 
duality. Our results demonstrate that the dual conformal symmetry persists at the quantum 
level up to an anomaly which is identical to that of a Wilson loop. We interpret this as strong 
evidence for the existence of a dual Wilson loop which remains to be constructed. 

We list a number of open questions for future work. A first one concerns the status of the 
dual conformal symmetry away from the origin of moduli space, e.g. in the Higgsed theory. 
As demonstrated in section (5), to lowest order in the masses (logarithmic accuracy), the 
Higgsed theory enjoys an exact dual conformal symmetry under which the masses transform 
in a nontrivial way. It is not clear whether this symmetry extends all the way into the moduli 
space; for one thing, the origin of the symmetry is mysterious and the original string theory 
argument in [28] does not apply in ABJM due to difficulties with the T-duality. As discussed 
in the main text, a key step here would be to settle this question for the tree amplitudes. 

We note that a 3-loop computation of the 4-point and 6-point amplitude in ABJM would 
probably be feasible with the same techniques, although a more sustained effort would be 
required. For instance, we expect only degree-3 transcendental functions in the result. Fur- 
thermore, the only divergences should be double-logarithms multiplying the 1-loop amplitude. 
Given the absence of overlapping divergences, the integration technology developed in section 
(6) might thus plausibly be sufficient. 

An interesting property of our remainder function Rq is that it does not vanish in collinear 
limits, contrary to the case in SYM4. In fact, it even diverges logarithmically in the 'simple' 
collinear limit (six point goes to five), this even though the five-point amplitude is zero. This 
does not violate any physical principle, since the ^4g ree and ^-gstiifted P re f ac t° rs do not have any 
pole in this limit. In the absence of a pole, there is no need for the amplitude to factorize into 
a product of lower-point amplitudes. In other words, the leading term in the collinear limit in 
ABJM is similar to subleading, power-suppressed terms in the collinear limits in D = 4. The 
factorization theory for these terms is more complicated, and in fact it has only been worked 
out recently in the dual Wilson loop language [47]. It would be very interesting to work out 

13 Note that the vanishing of 1-loop Wilson loops was obtained numerically in [42]. Given the subtle analytic 
properties discussed in section (3.3), it could be worthwhile to supplement this by an analytic computation. 
14 At least up to 5 3 (P)5 6 (Q) and a purely bosonic factor, akin to the Parke-Taylor denominator in SYM4. 
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the general structure of this limit, using field theory arguments, as this should place strong 
constraints on the amplitudes. In subsection (3.3), for instance, we have conjectured from 
analyticity of the scattering amplitudes that a certain discontinuity of the amplitude should 
vanish in the collinear limit, but we have no idea how this could be established. 

Also interesting are the 'double-collinear' limit (six point goes to four), or factorization 
limits (Pi23 goes to zero, but momenta pi,2,3 do not become collinear). Since it was not clear 
to the authors what kind of field theory predictions are available for these limits, we did not 
discussed them on our 6-point result. However, it is possible that this could shed further 
light on our result itself, for instance by giving a physical interpretation for the relative signs 
between different terms. These limits may also yield some interesting constraints on the 
higher-point amplitudes. 

Another interesting direction for future work concerns the rest of the Yangian algebra 
at loop level. As one easily sees from [3], the Yangian algebra in ABJM is generated by the 
bosonic dual conformal symmetry together with the (ordinary) superconformal symmetry. 
Since the former is presently conjectured to become anomaly-free to all loops after dividing 
by the BDS Ansatz, the crux is the superconformal anomaly. By analogy with SYM4, the 
properly understood symmetry at the quantum level should uniquely determine the ampli- 
tudes, providing for an efficient way to compute them. Our two-loop result (7.5) should thus 
provide an important data point to understand the quantum symmetries of ABJM, perhaps 
combining the 1-loop ABJM analysis in [14] with the all-loop SYM4 analysis in [48]. 
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A Identities 

In this appendix, we aim to prove a series of identities used in the text. First consider the 
following identity: 

Ci+Cj- _ e(6,l,2,3,4)(3-5)(5-l) 

C 2 +C* e(3,4,5,6,l)(6-2)(2-4)' 1 ' ] 

The strategy is to express the five-dimensional e symbol in terms of angle brackets. To do so, 
we start from the definition of the e-symbol as a determinant and use the manifest translation 
invariance of the formula to set %2 = 0. In doing so, we must remember to normalize the 
determinant such that e(i, j, k, I, m) 2 agrees with the Gram determinant formula (3.16), since 
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this is the convention used in the main text; this requires an extra factor of 2i\/2. Thus 
e(6, 1, 2, 3, 4) := 2»>/2 det(j/6, y u y 2 , y 3 , y 4 ) = 2i^2 det 



' -Pe-pi -pi p 2 P2+m^ 

1 1111 
^ -(61) 2 -(23) 2 y 



where the first three rows are real in Minkowski signature. This determinant can now be eval- 
uated in terms of three-dimensional ones, which in turn give two-brackets: det (pi,Pj,Pk) '■= 
\{ij){jk){ik) . This way we obtain 

e(6,l,2,3,4) = *V2(61)(12)(23)((31)(16) + (32)(26)). (A.2) 

Performing a similar computation for e(3,4, 5,6, 1) and using that (3-5) = — (34) 2 etc., we 
thus find 

e(6,l,2,3,4)(3-5)(5-l) = (12)(34)(56) / (31)(16) + (32) (26) \ 
e(3,4,5,6,l)(6-2)(2-4) (23) (45) (61) V (34) (46) + (35) (56) ) ' 1 ' 

Using momentum conservation, the parenthesis can be shown to equal — 1 , proving the desired 
formula using (3.10). 

Another remarkable algebraic identity is 

x\A - x\A = ((34) (46) + (35)(56)) 2 (A.4) 

which one might call a Dirac matrix trace identity, and follows from squaring eq. (A.2) and 
using that the square should give the Gram determinant. Using this identity we have that 

(12) (34) (56)6(3,4,5, 6,1) _ ^(12) (45) ((34) (46) + (35) (56)) 



2 



(l-4)(3-6)(l-5)(3-5)V«i(l-«i) ^(l2) 2 (45) 2 ((34)(46) + (35)(56)) 
which was used around eq. (7.2). 

B ABJM theory on the Higgs branch 

The action of ABJM takes the form (see for instance [49] for an explicit component form): 

£ = ^(£kin + A + A>). (B.1) 

To describe the spectrum of the theory on the Higgs branch, we begin by describing the 
fermion mass matrix. The interactions of the fermions can be written, following [49] but 
as can also be verified directly by comparing against various components of the four-point 
amplitude (2.7), 

U = (TrtVuV^^D] - Tr[^ D c ij B }) (26$6% - 6%dg) 

+e ABC D^ A ^ B (t> C ^ D ] + e ABCD Ti [iPaMc4>d] ■ (B.2) 
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As already mentioned in the main text, the moduli space (C 4 /Zfc) Ar of this theory is charac- 
terized by diagonal vacuum expectation values for the scalar fields. Let us denote the fields 
above (below) the diagonal with a plus (minus) superscript, so that (V^)^ = V' j4=F . As one 
can easily see from the action (B.2), the diagonal fermions remain massless while tpg and 
^pB+ m - x each other. Upon inserting a diagonal vev for the scalars, the mass term thus 
takes the form (V>^~, tp^Mttyg, ^ B+ ) T where My is the 8x8 Hermitian matrix 



Mt 



2(xx - yy)% - 6% (x-x - y-y) 2e A BCDX C y D 



\ 2e ABCD y c x D 2(yy - xx)^ - 5^(y-y - x-x) J ' 

In this appendix, (x,y) A := (vi,Vj) A will denote the diagonal vevs coupled to off-diagonal 
components under consideration. As one can verify M 2 = t$m 2 with 

m 2 = (x-x + y-y) 2 — 4x-yy-x, 

showing that all 8 off-diagonal fermions acquire the same mass. 
The scalar potential was described in detail in ref. [1], 

£ 6 = Tr&% A <p B c] <p c B } - ^Tr[<^ [A <^ B c c] ] 

where here the square bracket means antisymmetrization in the indices. Again one can see 
that the diagonal fluctuations remain massless while off-diagonal ones 5(f> A+ and 8(t>\ mix 
with each other. It follows that the mass term takes the form (S<pJ, 5(j) A -)M 2 (8(i) B+ , <5<^) T , 
and a computation gives the 8x8 Hermitian matrix as 

M 2 = m 2 (Is - P 8 ) (B.3) 



where = _ T _ T 1 • __ _ I • I /m 2 is an orthogonal pro- 

y-y 1 -x 1 J y -2y-xy x-x + y-y J \^y -x J 

jector onto the two would-be Goldstone bosons (5(p + ,5(f> + ) ~ (x,—y) and 

(y,—x). We conclude that six of the eight scalars acquire the same mass squared as the 

fermions, while the remaining two acquire no mass, although they are soon to be "eaten" by 

the gauge fields through the Higgs mechanism. 

Finally, we consider the gauge fields, which acquire mass terms through the scalar kinetic 

term ~ TrD^D^fi. First we discuss the off-diagonal components. As one can see again the 

fields Af 2 from the two gauge groups mix with each other, so the mass term is characterized 

by a 2 x 2 Hermitian matrix (A^[ , A 2 )M g (Af , A 2 ). However, in this case the kinetic term is 

also characterized by a nontrivial matrix d(A~^ , A 2 ) A K g (Af , A 2 ). These two matrices are 




and ,U 7 f • - 2X ^ 
\ —2y-x x-x + y-y 



Fortunately, to obtain the propagator it is not necessary to diagonalize these two matrices 
simultaneously — as pointed out in [50] this may not even be possible in general. In the 
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present case, one can verify that (K g M g ) 2 = m 2 l2 and this suffices in order to write down 
the propagator in a simple way. To see this, let us first add a gauge-fixing term to the action 
{d^A^ — ^v-(6<j)~))Kg(d l j_A ,1+ — £f-(£</> + ))/£, designed to remove the mixing between the 
gauge bosons and the scalar fields, where £ is some arbitrary scale. Then the two unphysical 
scalars acquire masses squared ~ £m, and a short computation gives the gluon propagator as 



' ' A ° )(p> U + J >0t i> + (W + ?? p' + (HW ' 

In particular, this formula shows that there are no singularities at zero momentum provided 
m 2 7^ 0, as required in the main text. 

Finally, we discuss the diagonal gauge fields. Since only the combination [A\ — A2) 
receives a mass term in this case, we have that M g oc (1, —1) <X> (1, — 1) T which is effectively 
nilpotent: (K g M g ) 2 = 0. As the above propagator shows, even though the mass matrix is 
nonzero, no massive states appear in the spectrum (as required by supersymmetry) . This 
situation has been discussed in detail in [50]. 15 

C Integrals using the mass regularization 

Here we summarize the results obtained in section (6) for the integrals defined in eqs. (4.5), 
multiplied by 16-7T 2 , evaluated using a small internal mass to regulate infrared divergences as 
defined in section (5). 



I (1) - + — + -(l + logli 3 Jlog 2 2 --lOg 2 2 
1 z x 13 x 15 ^ x 13 x 15 

1, 



„2 iJl \ t ; M „2 / 2 



-Li 2 (l - x 13 /x 36 ) - Li 2 (l - x 15 /x 25 ) + -Li 2 (l - l/u 3 ), 



r2mh/ n _ 9 ^ 1 1 2 £/fmgl5^35 r-,, ,.2 / 2 \ 2 2 

J even I 1 J ~ 2 ^ 7 io S 2 4 L1 2^ X 25/ X 15J - X 25 / £ 35 J, 

T 13 X 25 



rl3;13 _ _ , fMm 
x 13 

7 13;35 _ 1 _ 1 , 4 ^IR. a; 15 
z x 13 x 35 

^ m rf h (i) = o. (C.l) 

In addition, we have the following two absolutely-convergent integrals: 

I cri t tcr (1)+/ 13;46 = _ 1 ^ (J _ ^ _ 1 ^ (1 _ ^ _ 1 ^ ^ bg ^ _ (arccos( ^ T)) 2 + _^ (C2) 

15 Note that in refs. [50] it was further shown that the field {A\ — A%) can be integrated out in a systematic 
expansion in 1/m, yielding a Yang-Mills term kF 2 /m for the remaining gauge field plus other terms. But 
since for us m ~ \vi\ 2 is an infrared scale, not an ultraviolet scale, such a (in any case not strictly necessary) 
procedure would be inappropriate in our context. 
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and (see eq. (6.18) for the definition) 



rodd 

^boxitri 



1) 



log ^ arccos( 



>ul) 



(C.3) 



20x1(1 - ul) 

For completeness, we find for the four-point double-box using the same regularization: 

_2 



r l,2,3;3,4,l 
l 2box 



[e(a, 1, 2, 3, *)e(b, 3, 4, 1, *)] = ^ + log ^ 



log log (C.4) 



,r 



13 



21 



D Integrals using dimensional regularization 

In this appendix, we present the integrated result of infrared divergent integrals using di- 
mensional regularization. Here all integrals are again multiplied by 16-k 2 . After obtaining 
the integrals in terms of Feynman parameters, we integrate by converting the integrand into 
Mellin-Barnes representation, and implement the Mathemtica package MB.m [51] to obtain 
the result up to 0(e). The result is expressed in terms of zero-, one- and two-dimensional 
integrals in Mellin space. The one and two-dimensional integrals are analytically evaluated 
by performing sum over residues. That such sum can be carried out analytically, is simply 
due to the fact that the two-loop amplitude should be of transcendental two functions. 
The two mass hard integral gives: 



/ 2mh m 



e(o ) l,2,3,*)e(6,3,5,l*) 



(a ■ l)(o • 2)(a ■ 3) (a • b)(b ■ 3) (6 • 5) (6 • 1) 



-7 2e 



(4vr)- 2e 
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-2<En 
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log 2 (x( 3 /x( 5 ) - log'(x( 3 /xi 5 ) + 21og 2 (x 2 3 /* 2 5 ) 



+21og 2 (x 2 5 /xy -210^5/45) + 3 log 2 (^5/^5) 



2/™2 



,2/2 



+vrarcsin(A/x| 5 /x| 5 ) - arcsin (Jx^/xfe) + ir axcsm(J xfo/afa) - arcsin (Jx% 5 /x% 5 ) 



-Li 2 (l- x 2 15 /x% 5 ) 



1 . 2 . 2 . llvr 2 
-Li 2 {l-x 35 /x 25 ) + 



log 2 2 



2 - 2vr 2 * a 



(D.l) 
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where a « 0.3594267177020808. The crab integral gives: 



I 4 crab (l) 



e(a,l,2,3,*)e(6,5,6, 1,*) 



-7 2e 



(4vr)- 2e 
1 



>• l)(a-2)(o-3)(a-6)(6-5)(6-6)(fe-l 

(4»)- 2e + (^i 5 )' 2£ + O^r 26 - (4,r 2e - (■'■ 

8e 2 



36 ) _|_ ( x 35) 



2c 



(4»)- 2e 

4e 



(*? 5 ) 



-2eH 



!°g 2 (^13/^15) - V^ia/^l) - log 2 (^ 2 5 /^36) - log^s/slts) - log 2 (x 2 3 /a;^) 
W^le/aas) + log 2 (xl 6 /x 2 15 ) + log 2 (x^ 5 /xl 5 ) + log 2 (x 2 6 /x§ 5 ) - log 2 (x| 5 /xl 6 ) 



7rarcsin(Wxf 5 /a;| 5 ) + arcsin (\/x{ 5 /x% 5 ) - arcsm(Wxf 3 /x§ 6 ) + arcsin (Jxf 3 /x 



+ 7 



log 2 (u 3 ) + L*2(l - ^15/^25) + L *2(l - xl 3 /xl 6 ) + Li 2 (1 - u 3 ) 



TT Z 1 

+ 48 + 2 
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Notice the presence of arcsin functions with non-conformal cross-ratios as arguments. Such 
function did not appear in the mass regulated result and marks a stark distinction between 
the two regularizations. The ArcSin functions always come in the combination 

arcsin(y / m) — arcsin 2 (y/m)/TT . 



This particular combination is necessary for the integral to remain real. 
For completeness, we list the double triangle result: 



I. 



1,3;3,1 
2tri 



I, 



1,3;5,1 
2tri 



(1-3) 2 

(o-l)(o-3)(o-6)(6-l)(6-3) 
(l-3)(l-5) 



e 'x 



13 



-2c 



4tt 



13 •''IS 



4ttx 2 5 



1 

Ye 
-2e 



1 1 



(a-l)(o-3)(a-6)(6-5)(6-l) 

While the integrated results appears to be regularization scheme dependent (compare 
with eqs. (C.l)), when combined into amplitudes they give identical result up to additive 
constants. In particular, the arcsin functions completely cancel. Considering the sum of 
infrared divergent integrals in eq. (7.1) one obtains: 
6 



E ^n(*) + / Crab «+S 2;J+2 '* 
i=l 



6 

£ 

i=i 



-2c 



•in 



,-2e 



{2ef 



log 



l 2tri 



V ii+2 1 „ X i+li+3 



log- 



"M + 3 



Hi+3 



. 1 , 2 x m+3 

+ 7 log 



:Lio 1 



•71 2 ™ 2 

■ f/ ii+4 x i+li+3 

/y 2 , yi 2 

J 'ii+3 J 'i+li+4 

,31 



,23 



7T 



12 * a) 



5BS 6 (e -»• 2e) -vr 2 



12 * a) 



(D.3) 



where BDSq is the one-loop six-point MHV amplitude of = 4 sYM [30] with e replaced 
by 2e, reflecting the two-loop nature of the result. Thus the dimensionally regulated infrared 
divergent integrals combine to give the BDS answer, just as in the mass regulated result. 



- 41 - 



References 

[1] O. Aharony, O. Bergman, D. L. Jafferis and J. Maldacena, gravity duals," JHEP 0810, 091 
(2008) [arXiv:0806.1218 [hep-th]]. 

[2] T. Bargheer, F. Loebbert and C. Meneghelli, Phys. Rev. D 82, 045016 (2010) [arXiv:1003.6120 
[hep-th]]. 

[3] Y. -t. Huang and A. E. Lipstein, JHEP 1011, 076 (2010) [arXiv: 1008.0041 [hep-th]]. 

[4] J. M. Drummond, J. Hcnn, G. P. Korchcmsky and E. Sokatchev, Nucl. Phys. B 828, 317 (2010) 
[arXiv:0807.1095 [hep-th]]; 

[5] J. M. Drummond, J. Henn, V. A. Smirnov and E. Sokatchev, JHEP 0701, 064 (2007) 
[hep-th/0607160]; 

J. M. Drummond, J. M. Hcnn and J. Plefka, JHEP 0905, 046 (2009) [arXiv:0902.2987 [hep-th]]. 

[6] D. Gang, Y. -t. Huang, E. Koh, S. Lee and A. E. Lipstein, JHEP 1103, 116 (2011) 
[arXiv: 1012.5032 [hep-th]]. 

[7] N. Bcrkovits and J. Maldacena, JHEP 0809, 062 (2008) [arXiv:0807.3196 [hep-th]]; N. Beisert, 
R. Ricci, A. A. Tseytlin and M. Wolf, Phys. Rev. D 78, 126004 (2008) [arXiv:0807.3228 
[hep-th]]. 

[8] I. Adam, A. Dckcl, Y. Oz, JHEP 0904, 120 (2009). [arXiv:0902.3805 [hep-th]], JHEP 1010, 110 
(2010) [arXiv: 1008.0649 [hep-th]]; 

P. A. Grassi, D. Sorokin, L. Wulff, JHEP 0908, 060 (2009). [arXiv:0903.5407 [hep-th]]. 
I. Bakhmatov, Nucl. Phys. B 847, 38 (2011) [arXiv: 101 1.0985 [hep-th]]. 

[9] N. Arkani-Hamed, F. Cachazo, C. Cheung and J. Kaplan, JHEP 1003, 020 (2010) 
[arXiv:0907.5418 [hep-th]]. 

[10] S. Lee, Phys. Rev. Lett. 105, 151603 (2010) [arXiv: 1007.4772 [hep-th]]. 

[11] A. Agarwal, N. Beisert and T. McLoughlin, JHEP 0906, 045 (2009) [arXiv:0812.3367 [hep-th]]. 

[12] Y. -t. Huang, "Recent Advances in Scattering Amplitude" INI Cambridge, 
http://www.newton.ac.uk/programmes/BSM/seminars/040409001.html 

[13] M. S. Bianchi, M. Leoni, A. Mauri, S. Penati and A. Santambrogio, arXiv: 1204.4407 [hep-th]. 

[14] T. Bargheer, N. Beisert, F. Loebbert and T. McLoughlin, arXiv: 1204.4406 [hep-th]. 

[15] A. Brandhuber, C. Wen and G. Travaglini, arXiv:1205.6705 [hep-th]. 

[16] W. -M. Chen and Y. -t. Huang, JHEP 1111, 057 (2011) [arXiv:1107.2710 [hep-th]]. 

[17] M. S. Bianchi, M. Leoni, A. Mauri, S. Penati and A. Santambrogio, JHEP 1201, 056 (2012) 
[arXiv:1107.3139 [hep-th]]. 

[18] A. Brandhuber, G. Travaglini and C. Wen, arXiv:1207.6908 [hep-th]. 

[19] Z. Bern, J. J. M. Carrasco and H. Johansson, Phys. Rev. D 78, 085011 (2008) [arXiv:0805.3993 
[hep-ph]]. 

[20] Z. Bern, J. J. M. Carrasco and H. Johansson, Phys. Rev. Lett. 105, 061602 (2010) 
[arXiv: 1004.0476 [hep-th]]. 



- 42 - 



Z. Bern, T. Dcnncn, Y. -t. Huang and M. Kiermaier, Phys. Rev. D 82, 065003 (2010) 
[arXiv: 1004.0693 [hep-th]]. 

J. Bagger and N. Lambert, Phys. Rev. D 77, 065008 (2008) [arXiv:0711.0955 [hep-th]]; 

Phys. Rev. D 79, 025002 (2009) [arXiv:0807.0163 [hep-th]]; 

A. Gustavsson, Nucl. Phys. B 811, 66 (2009) [arXiv:0709.1260 [hep-th]]. 

T. Bargheer, S. He and T. McLoughlin, arXiv: 1203.0562 [hep-th]. 

N. E. J. Bjerrum-Bohr, P. H. Damgaard and P. Vanhove, Phys. Rev. Lett. 103, 161602 (2009) 
[arXiv:0907.1425 [hep-th]]. 

M. S. Bianchi, M. Lconi and S. Penati, JHEP 1204, 045 (2012) [arXiv:1112.3649 [hep-th]]. 

J. M. Drummond, J. Henn, G. P. Korchemsky and E. Sokatchev, Nucl. Phys. B 795, 52 (2008) 
[arXiv: 0709. 2368 [hep-th]]. 

A. Brandhuber, P. Heslop and G. Travaglini, JHEP 0908, 095 (2009) [arXiv:0905.4377 
[hep-th]], JHEP 0910, 063 (2009) [arXiv:0906.3552 [hep-th]]. 

L. F. Alday, J. M. Henn, J. Plcfka and T. Schuster, JHEP 1001, 077 (2010) [arXiv:0908.0684 
[hep-th]]. 

Z. Bern, L. J. Dixon and V. A. Smirnov, Phys. Rev. D 72, 085001 (2005) 
[arXiv:hcp-th/0505205]. 

Z. Bern, L. J. Dixon, D. C. Dunbar and D. A. Kosower, Nucl. Phys. B 425, 217 (1994) 
[hcp-ph/9403226]. 

P. A. M. Dirac, Annals Math. 37, 429 (1936); 

G. Mack and A. Salam, Annals Phys. 53, 174 (1969); 

S. L. Adler, Phys. Rev. D 6, 3445 (1972) [Erratum-ibid. D 7, 3821 (1973)]; 
R. Marnelius and B. E. W. Nilsson, Phys. Rev. D 22, 830 (1980). 

W. Sicgel, arXiv:1204.5679 [hep-th]. 

W. L. van Necrvcn and J. A. M. Vcrmascrcn, Phys. Lett. B 137, 241 (1984). 

R. Britto, F. Cachazo and B. Feng, Nucl. Phys. B 715, 499 (2005) [hep-th/0412308]. 

Y. -t. Huang and S. Lee, arXiv:1207.4851 [hep-th]. 

Z. Bern, J. S. Rozowsky and B. Yan, Phys. Lett. B 401, 273 (1997) [hcp-ph/9702424]. 

J. L. Bourjaily, A. DiRe, A. Shaikh, M. Spradlin and A. Volovich, JHEP 1203, 032 (2012) 
[arXiv:1112.6432 [hep-th]]; 

J. Golden and M. Spradlin, JHEP 1205, 027 (2012) [arXiv:1203.1915 [hep-th]]. 

N. Craig, H. Elvang, M. Kiermaier and T. Slatyer, JHEP 1112, 097 (2011) [arXiv:1104.2050 
[hep-th]]; M. Kiermaier, arXiv: 1105.5385 [hep-th]. 

S. Caron-Huot and D. O'Connell, JHEP 1108, 014 (2011) [arXiv:1010.5487 [hep-th]]. 

T. Denncn and Y. -t. Huang, JHEP 1101, 140 (2011) [arXiv:1010.5874 [hep-th]]. 

W. Chen, G. W. Semenoff and Y. S. Wu, Phys. Rev. D 46, 5521 (1992) [arXiv:hep-th/9209005]. 

J. M. Henn, J. Plclka, K. Wiegandt, JHEP 1008, 032 (2010) [arXiv: 1004.0226 [hep-th]]. 

K. Wiegandt, Phys. Rev. D 84, 126015 (2011) [arXiv:1110.1373 [hep-th]]. 



-43- 



[44] M. F. Paulos, M. Spradlin and A. Volovich, JHEP 1208, 072 (2012) [arXiv:1203.6362 [hep-th]]. 
[45] S. Caron-Huot, K. J. Larsen and K. J. Larsen, arXiv:1205.0801 [hep-ph]. 

[46] L. J. Mason and D. Skinner, JHEP 1012, 018 (2010) [arXiv: 1009.2225 [hep-th]]; S. Caron-Huot, 
JHEP 1107, 058 (2011) [arXiv:1010.1167 [hep-th]]; see also, B. Eden, P. Heslop, 
G. P. Korchcmsky and E. Sokatchev, arXiv:1103.3714 [hep-th] and arXiv: 1103.4353 [hep-th]. 

[47] L. F. Alday, D. Gaiotto, J. Maldaccna, A. Sever and P. Vieira, JHEP 1104, 088 (2011) 
[arXiv: 1006.2788 [hep-th]]. 

[48] S. Caron-Huot and S. He, JHEP 1207, 174 (2012) [arXiv:1112.1060 [hep-th]]. 

[49] J. A. Minahan and K. Zarcmbo, JHEP 0809, 040 (2008) [arXiv:0806.3951 [hep-th]]. 

[50] S. Mukhi, JHEP 1112, 083 (2011) [arXiv:1110.3048 [hep-th]]; S. Mukhi and C. Papageorgakis, 
JHEP 0805, 085 (2008) [arXiv:0803.3218 [hep-th]]. 

[51] M. Czakon, Comput. Phys. Commun. 175, 559 (2006) [hep-ph/0511200]. 



- 44 - 



